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CHAPTER I. 


GENERAL THEORY 


1. — Thermodynamical Theory of Irreversible Processes. 


The thermodynamical theory of irreversible processes (7?) is a relatively 
new branch of science which enables us to treat real processes, from a macro- 
scopic point of view, in a systematic way. The theory is based on two fund- 
amental concepts, the entropy production and the Onsager relations (*), which 
are briefly introduced in this section. 

The entropy production, which is a result of the progress of irreversible 
processes, can be calculated with the help of the fundamental laws of macro- 
scopic physics: the law of conservation of mass, the equations of motion 
(momentum law), the first law of thermodynamics (energy conservation law) 
and the second law of thermodynamics (Gibbs’ relation). From the fund- 
amental laws a balance equation for the entropy can be derived: 


(1.1) ds,/ ot = — div (J, + sv) +0, 


which shows that the entropy s, (per unit volume) changes as a result of an 
entropy flow J, and an entropy production o. The divergence of the entropy 
flow gives the entropy supplied by exchange of energy or matter with the 
surroundings (it may be positive or negative), while o arises from the action 
of irreversible processes inside the system. Calculations show that the en- 
tropy production has an expression of the following form: 


(1.2) (oy di J; ne = 0 ; 


where we have a sum of products of quantities J; called «fluxes » (e.g., heat 
flow, diffusion flow, electric current, chemical reaction rate) and corresponding 
quantities A, called « forces» or «affinities» (e.g., temperature gradient, gra- 
dient of chemical potential, electric field strength, chemical affinity). This 
expression is a measure of the irreversibility of the processes: it vanishes when 
the system undergoes only reversible changes, in all other cases it is positive. 

In the most general case every flow is caused by contribution of all forces 
(and inversely). So in the linear approximation the irreversible phenomena 


‘an be expressed by the so-called « phenomenological equations », which have 
the following form: 


(1.3) J; = De het . 


() S. R. DE Groot: Thermodynamics of Irreversible Processes ( 


Amsterdam and 
New York, 1951). 


(2) I. PRIGOGINE: Étude thermodynamique des phénomènes irréversibles (Liege, 1947). 
(5) L. OxsaGer: Phys. Rev., 37, 405 (1931): 38, 2265 (1931). 4 ; 
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Examples of these laws are, for instance, Fourier’s law for the heat conduction, 
the law which describes the Dufour effect, Fick’s law for the diffusion, the 
thermal diffusion (Soret-effect) law, Ohm’s law for the electrical conduction, 
the laws which describe the thermoelectric effects and the laws which describe 
the galvanomagnetic and thermomagnetic effects. 

The coefficients D,, are called « phenomenological coefficients », and ONSA- 
GER (8) has established the following «reciprocal relations » which hold bet- 
ween them: 


(1.4) Lin = Lx, 


expressing a connection between reciprocal phenomena which arise from mutual 
interference of simultaneously occurring irreversible processes (*). Many well- 
known physical relations, which could not be explained by thermostatie con- 
siderations, or by considerations on the physical symmetry of the system, 
are a consequence of the « Onsager relations ». It is allowed to say that the 
Onsager relations are indeed a new and very fundamental law of macroscopic 
thermodynamics. ONSAGER has shown that they are the macroscopic counter- 
part of the property of time reversal invariance of the laws which govern the 
motion of the individual particles forming the system. The demonstration 
of ONSAGER is based on the general ideas of statistical mechanics, and it will 
be explained in a modified form in the next section. 


9. — De Groot and Mazur’s Theory. 


2:1. Casimir’s objection and method. — As CASIMIR (‘) pointed out, ONSA- 
GER’s proof is strictly valid only for scalar phenomena. As a matter of fact, 
ONSAGER assumes that the irreversible fluxes can be considered time deri- 
vatives of thermodynamic state variables. This is correct for scalar processes 
(such as chemical reactions and relaxation phenomena) but not for vectorial 
processes, (such as heat conduction, diffusion and electrical conduction), and 
tensorial processes (such as viscous flow), and therefore an extension of the 
theory is necessary. 

CASIMIR has proposed a method, which he has applied to the case of heat 
conduction in anisotropic crystals without magnetic field: he reformulates the 
problem under investigation in such a way that it is possible to apply the 
original Onsager formalism. To achieve this aim it is necessary to introduce 
appropriate auxiliary phenomenological equations containing an infinite number 


(*) In case there is a magnetic field B, the phenomenological coefficients can 
be functions of B and then the relation (1.4) must be replaced by L(B) = L,,(—B). 
(4) H. B. G. Casimir: Rev. Mod. Phys. 17, 343 (1945). 
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of coefficients, which are only indirectly related to the ordinary macroscopic 
laws and coefficients. Mazur and DE Groor (5) have also considered systems 
with magnetic fields, and derived reciprocal relations for the heat and electrical 


conduction tensors. 


9:2. Basis of the method of De Groot and Mazur. — To solve the problem 
of finding reciprocal relations for vectorial and tensorial processes, DE GROOT 
and MAZUR (*) have developed a different method, which is in a certain sense 
the opposite of the first: now the phenomenological problem is not rewritten 
in a form appropriate for the application of the original Onsager theory, but 
the formalism is generalized in such a way that it can be applied also for fluxes 
which describe vectorial and tensorial irreversible phenomena, or, in other 
words, for fluxes which are not necessarily time derivatives of state variables. 
In this theory the ordinary macroscopic laws can be used as phenomenological 
equations, and thus the complications inherent in the approach of the first 
method are avoided. With the help of a generalization of the fluctuation theory, 
(cf. also (°"*)), reciprocal relations for vectorial and tensorial processes can be 
derived directly from the property of microscopic reversibility. 


2°3. Fluctuation theory. — We consider an energetically insulated system. 
The thermodynamic state of the system is described by variables A,, Ag,..., Am, 
which are even functions of the particle velocities, and variables B,n413 Bris, 

, B,, which are odd functions of the particle velocities (the B-type variables, 
or ee divergences, may be time derivatives of A-type variables: cf. also (*)). 
All variables are continuous functions of space and time coordinates. For 
convenience the system will be divided into a number of cells of volume V", 
in which the thermodynamic variables may be considered as uniform (here 
4 numbers the cells). The deviations of the state variables from their equi- 
librium values are denoted by 


(1.5) aff = AM — (AM), (i= 1,2, ..., m) 
(1.6) (ee iy D) (k= ml, eee en) 


Since the entropy has a maximum in the equilibrium state, the deviation of 
the entropy of the system from its equilibrium value is given, as a first appro- 
ximation, by the quadratic form 


(127) AS=—-1 Da Ae ox grote! + Sr hit? Be a) 


where gf and hi’ are positive definite forms. 


(5) P. Mazur and S. R. pe Groot: Physica, 19, 961 (1953). 
(5) S. R. DE Groor and P. Mazur: Phys. Rev., 94, 218 (1954); 
Supplemento al Nuovo Cimento, 12, 5 (1954). 


(5:bis) J. VLIEGER and §. R. DE GROOT: Physica, 20, 372 (1954). 
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Since the entropy is an even function of the particle velocities no cross-terms 
between «- and f-variables appear in the expression (1.7). The probability 


IT 


distribution for the x‘ 


t 


and pi! is expressed by 


(1.8) TRARRE [Idata 
|... | exp[AS/E] Tux dowd pe 


where k is Boltzmann’s constant. 
The following linear combinations of parameters are introduced: 


(1.9) RIVE GAS ca —= VS qu, 
(1.10) a Vo AS ope > aD 


which will be designed as conjugate variables. 
The following averages can easily be found with the help of (I.8)-(1.10): 


(1.11) aeli(t) XY (t) == —- bi; On| V" 


(I.12) AU) 1 (6) = du Ov) V 4 , 


where 6 is the Kronecker symbol. 

Here and in the following, according to the fundamental notions of sta- 
tistical mechanics, the average may be interpreted either as an average over 
a microcanonical ensemble of systems, or as a time average for one single system. 

Average values of products of «- and B-type variables always vanish. 

Passing to the limit of continuous variables, the indices y and v of (1.11) 
and (1.12), which indicate the cells, can be replaced by r and r’, which denote 
space coordinates, whereas the last two factors of (1.11) and (1.12) combine 
into a Heaviside-Dirac 6 function. Consequently (I.11) and (1.12) become 


(1.13) ar, t)X;(r’, t) = — kd, d6(r—r’), 
(1.14) bi(r, tyr’, D = — kod,,6(r—r7’) . 


From (1.13) and (1.14) one can immediately derive the formulae 


(1.15) ar, HOQ(F')X (n't) = — bb, Q(r') Ô(r— r'}, 


(1.16) Br, HQ(r’)yi(r’,t) = — À Oy, Q(r') O(r — r'), 
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where Q(r) is a differential operator of the general form 


(1.17) Q(r) = Fg Ung sr) PTT [0A Oa, OM, « 


The coefficients a,,, are independent of the state variables f,(r,t) and «,(r, t); 
the Cartesian coordinates arc denoted by 2, 23 and 43. 

The formulae (1.15) and (I.16) serve as an extremely useful basis for the 
straightforward derivation of the reciprocal relations among phenomenological 
coefficients. 

Remark. It may be noted that with (1.9) and (1.10) the time derivative 
of (1.7) (entropy production per unit time) can be written as 


(1.18) AS = >, V" (D, axe Da Beye), 
or 
(1.19) ifs | {Dax (0) +d, Batr)yalr)} dr , 


in the limit of infinitely small cells. One sees that state variables and conjugate 
variables can be identified either from the expression for the deviation of the 
entropy, or from the expression for the entropy production (1.19). 


2°4. Microscopic reversibility. — The fact that, on the average, the future 
behaviour of a system is identical with its past behaviour, will now be expressed 
with the help of correlation functions. Since a magnetic field B is present, 
it is necessary to reverse its direction in every point in order to have the 
particles retracing their paths. Minus sign arises when an «-type variable 
is combined with a f-type variable (*"): 


(1.20) alr, tha, (r',t+-7){B, B= a, (r', t)x,(r,t+7){_ B, —B’, 
(1.21) adr, 4p (ry ET) {By BY = Br tar ET) By, Be 
(1.22) Bir, t)Pn(r”, t+ T) {B, B’; = Blt 0), ( (r, tt) {= B, iva B'}. 


From these three relations it follows that 


(1.23) ar, t)(0/0)x; (7,1) {B, B} = a; (r’, t)(0/ot)a,(r, t) {— B, — BY}, 
(1.24) afr, (0/00)f. Cr 0 {B, By = —B, (7, Dan) B,-—B%, 
(25) Bi(r, t)(0/et)B,.(r’, t) {B, BY = B(r!, t)(0/0t)Bi(r, t) (— B,— B. 


B and B' indicate the magnetic field Strenghts at the positions r and r' 
where the averages are performed. The time derivative in the last three re- 
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lations should be interpreted on a microscopic scale as the average decay of 
the fluctuations, which is described by difference quotients (*'); however, for 
all practical macroscopic purposes this quotient can be considered as a real 
derivative. 


2°5. Regression of fluctuations. — In order to find the reciprocal relations 
between the phenomenological coefficients, it is necessary to express the time 
derivatives of the state variables x, and /, as functions of the other state 
quantities X; and y, or of their derivatives with respect to space coordinates. 
(As a matter of fact we shall write instead of the state variables, their de- 
viations from the equilibrium value). That can be achieved with the use of 
the conservation laws and the entropy balance equations, and one obtains 
expressions of the kind 4 


(1.26) Get, {Ct OG D MOV, 
(1.27) Pe DEL DE VE 


The operators 2,, and Q,, turn out to be like the following: L,,: , L;,-grad, 
div (Lys), div (L;-grad), ete. The subscripts indicate here the fluxes and the 
forces, not the Cartesian components). 


2-6 Derivation of the Onsager relations. — With the preceding results, and 
the usual assumption that the average decay of fluctuations follows the pheno- 
menological macroscopic laws, (which is true for linear phenomena, so that 
we can take for the average decay of fluctuations, which appears in formulae 
(1.22)-(1.23), the expressions (1.26) and (I1.27)), the derivation of Onsager 
relations is straightforward. Introducing (1.26) and (1.27) into (I.23)-(I.25) 
one has 


(1.28) TREIA, (e, BY) yar YB B= 
r'){ >, gr, — B)xalr) + Da ie — B)yalr Sg 
(1.29) r) {D, fr, Bur) + 3, 200°, B )atr){B, B= 
BIS. 2 (r,— B) xx(r) + >, 22%, — B)y(r)}{ 
= fr) {Dry B)xe(r) + 2,9 sr, — B)y(n)}{— B, —B’}, 
(1.30) Br) {S, QE’, BY) xe(r') + 3, 8h (1', BY) B, BY = 


ce GDS OF. ir; Vede Dai sa — B) y,(r)}{—B, a B'}. 
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According to (1.15) and (1.16) these relations lead to 

ea bones (ri, Bor =r kop (r, Bor re), 
32) —kb,.25%(r', B'hô(r— r') =4+ 6, Dr, — B)d(r—r’), 


(eso ee 60,2 ir’, B O(n 1) kom QE (r, — B)d(r—r'). 
Elimination of the Kronecker 6’s and the 6-functions gives now the relations 
for the phenomenological coefficients, which are contained into the operators Q, 
either simply 


(1.34) Lix(B) = ee B) ’ 
or in a form containing differential operators, e.g. 


(1.35) Div L,,(B) = Div L,.(— B), 


according to the actual expression of £ for the various cases. 


Remark. Onsager’s demonstration was valid for the case J; = 0x,/0t and 
Xi= X;. Reciprocal relations are then immediately obtained: since Q,, = Ly, 
and from (1.3), (1.23) and (1.15), follows Lj, = Ly. 


CHAPTER IT. 


DERIVATION OF THE RECIPROCAL RELATIONS 


1. — Introduction. 


In this chapter the general theory outlined above (5) is applied to the case 
of mixtures of n charged components in an electromagnetic field (7). In such a 
general case the irreversible phenomena of heat conduction, diffusion (including 
electric conduction) and cross-effects, and of viscous flow, occur. 

First the macroscopic phenomenological theory of the case mentioned is 
developed: in § 2 the entropy production in local form is calculated from the 
fundamental laws; in § 3 the phenomenological equations are written down. 


(7) R. Frescut, 8. R. DE Groot, P. Mazur and J. VLIEGER: Physica, 20, 67 (1954). 
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Subsequently the fluctuation theory is developed and the reciprocal re- 
lations for the phenomenological coefficients are obtained (§§ 4-10). In § 4 
we derive an expression for the total entropy production in an energetically 
insulated system in terms of fluctuations of independent thermodynamic 
state variables; in § 5 fluctuation averages are calculated, and in the next 
section (§ 6) the microscopic time reversibility property is expressed by means 
of correlation functions. From the phenomenological equations, the consery- 
ation laws and the entropy balance, we derive in § 7 the expression for the 
regression of fluctuations of the state variables. The equation for the decay 
of the vector potential requires special attention. In $$ 8, 9 and 10 the reci- 
procal relations for the phenomenological coefficients are derived. It is found 
that one gets physical results from time reversibility relations, in which both 
even and odd variables oceur (?). 

The last part of the chapter ($ 11 and 12) shows how the description of 
simpler phenomena can be obtained by specializing the results of the more 
general case. In particular the equations which describe the galvanomagnetic 
and thermomagnetic phenomena are considered. 


2. — Entropy Production in Local Form. 


Let us consider a mixture of n charged non-reacting components in an 
electromagnetic field and subject to external conservative (non-electrical) 
forces. Magnetization and polarization are neglected. The entropy production 
can be calculated from the fundamental laws of macroscopic physics (7°) 
(cf. chap. I, § 2). 


The law of conservation of mass. — As chemical reactions are not taken 
into consideration, this law can be written as follows for the component k: 


(IL.1) o de,/dt = — div, , (= 152, Rev) 


where 0 is the total density, ¢ the mass fraction of component k (> c, = 1) 
k=1 


and J; = 0x(v,— v) the diffusion flow of k with respect to the centre of gravity 
motion (0; = ¢x@ is density of component k, vx is velocity of k and v is bary- 


centric velocity). From the definition of J, it follows tbat ÿ Je Oe 


k=1 


thermore we have 


(IL.2) o dv/dt = div v 


(8) P. MAZUR and I. PrIGOGINE: Acad. Loy. Belg., Classe des Sciences, Mém., 
98, fasc. 1 (1953). 
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Here v = o is the specific volume. In (II.1) and (IL.2) d/dt is the substantial 
time derivative with respect to the centre of mass motion: 


(11.3) d/dt = 0/0t + v-grad . 
With (11.3) the formulae (II.1) and (II.2) can also be written as 
(IL.4) dox/0t = — div 9, v— div Jr, 
(11.5) d0/ot = — div ov. 
The momentum law. 


(11.6) o dv/dt =— Div P'+ So + >Y oF - 


k—1 k=1 


F is the Lorentz electromagnetic force per unit mass, acting on compo- 
nent k: 


(II.7) FP =e{E+c*(v,\B)}, (ha ty Deca) 


where E is the electric field strength, B the magnetic induction, ex the charge 
per unit mass of component k. F, is an external (non-electrical) conservative 
force, also per unit mass of k: 


(IL.8) F,,=— gradw,; with 0w,.jo=0, (ie == Ls 2s) 


where w, is the potential energy per unit mass of k. P is the pressure tensor 
(IT.9) P=pò+I, 


with ò the unit tensor, p the hydrostatic pressure and II the viscous pressure 
tensor. The sign + means transposing of the Cartesian components of a tensor. 
Introducing the total charge density 


(IL.10) 06 = > oxen, 


k=1 
and the conduction current 


(11.11) i= x ext, 


k=1 


formula (II.6) can be written as 


(11-12) 9 dv/dt = — Div Pt + ge{E + (A B)} + (GAB) + JS QE. 
k=1 
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Sealar multiplication of (11.12) with v gives the balance of the kinetic energy 
per unit of mass i = }v?: 


(11.13) o dl/dt = — v- Div Pt + peu E — ci: (vA B) + > ow F, . 


k=1 
Using (IJ.3), (11.5) and (II.8), we can write (II.12) alternatively as 
(II.14)  dg/0t=— Div (gv + P*) — pel gradg +  (0A/0t — v\B)} + 


+eti\B)— > ox grad w;, 


k=1 
where g = ov is the momentum density; furthermore use has been made of 
(11.15) E =— grady—c dA/0, 

where œ is the scalar potential and A the vector potential. 


Energy law. — The energy of the electromagnetic field per unit of volume is 
(11.16) jp = #(E* + B?)- 
The balance equation of f, is Poynting’s law 
(11.17) of, /et = — div c(E\B)— oev:E— i-E, 


where c(E/\B) is the Poynting vector. With (II.2) and (II.3) formula (LEI) 
can also be written as 


(11.18) o df/dt = — div {e(E/ B) — hv}— cev:E— i-E, 


where f= f, v is the electromagnetic energy per unit mass. 


n 


For the potential energy w — > cio per unit mass, we have 
k=1 


n 


(11.19) dw/dt = — div (> wd) — > on. F 


k=1 k=1 


where (II.1), (11.3) and (II.8) have been applied. 
The energy equation can be written as 


(11.20) o de/dt = — div J, , 
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where e is the total energy per unit mass and J, the energy flow. These are 


given by 


(11.21) e=u+l+f+uw, 


sat I v:P ! c(E\B)--f,v + > wedi; 


k=1 


(11.22) J 


e 


with « the internal energy per unit of mass. Equation (II.22) defines the heat 
flow J,. [The balance equation of e, = eo can be written as 


(11.23) dele = div. J, tà — dio) 
with 


(II.24) J, gp = Woe Lv 4 ww A Jc P + cE ABY+ > we, 


k=1 


which follows from (II.22) and e, =u, +l, +f, + w, (the index v means per 
unit volume). In the formula (II.24) the electromagnetic (« anti »)-convection 
term —f,v has disappeared, in place of which one has the ordinary con- 
vection terms u,v, l,v and w,v]. 

The first law of thermodynamics (internal energy balance) is found by 
deducting (11.13), (IT. 18) and (11.19) from (11.20) (with (II.21) and (II.22)): 


(11.25) . o du/dt =— div J,+ i: {E+ c(v\ B)}+4+ > F,-J,— P: Grad v. 
k=1 
Gibb’s law. 
(11.26) T ds/dt = du/dt + p dv/dt — Yu, de,/dt , 


k=1 


where T is the temperature, s the entropy per unit mass and w, the chemical 
potential per unit mass. 


Entropy balance. — Introducing (11.1), (11.2) and (11.25) and using also 
(IT.9) and (11.11), equation (11.26) can be written as a balance equation, 


(11,27) ods/dt= — div J, +o, 
where J, is the entropy flow 


n 


(11.28) = (LS pid IT; 


k=1 
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and where the entropy production o follows from 


(29) fos — J,‘ (grad T)/T+ > Ji [E + ev NB)} + F,— 


k=1 


— T grad (wx/T)] — IL: Grad v 0. 


The three terms on the right-hand side of (II.29) are the results of entropy 
production due to heat conduction, diffusion and viscous flow respectively. 

Hereunder we shall use another form of To, which follows from (11.3), 
(II, 15) and (II.28): 


(11.30) To=— J,-grad T — > Ji {grad ff, + ex 1(0A/0t — v \ B)} — 
k=1 


= Grad. te One 


where the quantity 


(11.31) Ur = Ur À EP | 10% (= 15 2;0,%) 


has been introduced. 
Finally from (II.3), (11.5) and (II.27) one finds 


(II.32) NOR J.) +0, 
where s,= os is the entropy per unit of volume. 


Remark. From (11.30) (or (II.29)) one can see that, when the barycentric 
velocity vanishes, the magnetic part of the Lorentz force does not give any 
explicit contribution to the entropy production. It can however be shown 
(in an analogous way as it has been done by HoovMAN and others (°) for the 
case of the Coriolis forces) that an implicit dependence may exist. 


3. — Phenomenological Equations. 


In (11.30) To is seen to be a sum of the products of two sets of variables, 
viz. the first set J,, J, and II, which are called «fluxes », and the second set 
1(04/0t— v /\B)} and — Grad v, which are called 


— grad T, — {grad fix + exe 
With these fluxes and forces we can 


«forces» or «affinities » (*) (cf. (1.2). 


(2) G. J. Hoyman, H. Hozran Jr., P. Mazur and $. R. pe Groot: Physica, 19, 
1095 (1953). 

(*) The parameter 7 which multiplies 
aken into account in the expression of the forces. 


c can be included in the phenomenological 


coefficients, hence it is not t 


2 — Supplemento al Nuovo Cimento. 
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establish the « phenomenological equations » as the linear relations (cf. (1.3)) 


(11.33) J, =—L,-grad T— Y Ly -{grad fi, + ec (0A/ot— v /\B)}, 
1=1 

(11.34) J, =—L,,:grad T— > Lx {grad & + e,c-1(0A/ot— v \ B)}, 
1=1 


3 
(11.35)  ;=— Y Lismn(Grad V)mn, (1,7 = 1, 2, 3) 
m,m=1 

(k and J denote chemical components, î, j, m and n Cartesian components). 
In these formulae cross-effects between fluxes and forces of different tensorial 
character are not taken into consideration, although such effects might exist 
in anisotropic media. The present formalism could of course be easily extended 
to such a case. (In isotropic media these cross-effects do not exist, according 
to the so-called Curie-principle). 

The tensors L,, L,, and L, (k,{=1,2,...,%) are not independent, be- 
cause the relation 


(11.36) SJ 0 
k=. 


holds. Equations containing independent coefficients, however, are easily 
obtained when the dependent flow J, is eliminated from the expression of 
the entropy production. From (11.30) and (11.36) one has 


n-1 
(LEST) To =— J,-grad T— > Ji: {grad (üi— fn) + (ex — €n)c (0A) ot — 
k=1 


— v/\B)}—I:Grad v, 


hence the phenomenological equations 


n=1 . 
(11.38) J,=— L,,:grad T— Y L, {grad (f,— fin) — (e, — en)c-*(0A/0t— v \ B)}, 
x 1=1 
n--1 
(11.39) J,=— L,,-grad T— Y L,,-{grad (@,— %,)— (e, é,)e(CA/ot—v \B)}. 
1=1 


(eli co net) 


By using (II.36) and by comparing the new sets of equations with equations 
(11.33) and (II.34) one can easily prove that the new coefficients coincide with 
those coefficients of the first set, which have the same subscripts, and it can be 
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shewn that the following relations exist: 


(II. 40) DD UD = 0 2D Ue 0 > Le<0. (61,240) 
k=1 k=1 k=1 k=1 
We now consider the phenomenological equations (11.35). We shall take 
the case usually considered of a symmetric tensor II. Then in (11.30) only 
the symmetric part of the tensor Grad v is left. Therefore, instead of (11.35), 
we have 


3 
(II.41) ie. oe > L,;5, (mn (Grad V)cmn) , 

m,n= 1 
where the notations (ij) and (mn) indicate the symmetrical character in the 
indices between brackets. Consequently we are left with 36 phenomenological 
coefficients DL, ;;)o.n instead of 81 coefficients L;,,, from (11.35). 


4. — Entropy Production in Terms of Fluctuations. 


In order to apply the fluctuation theory one should know the fluctuation 
of the entropy (1.7) of the whole energetically insulated system, in terms of 
the fluctuations of the parameters (1.5), (1.6) which determine the thermodyna- 
mical state of the system, and of the variables conjugated to these para- 
meters (1.9), (1.10). However one can see that state variables and conjugate 
quantities can as well be found from the expression (I.19) for the entropy 
production per unit time, which is generally used instead of the fluctuation 
of the entropy (*). 

The change of the entropy per unit time dS/dt of the whole adiabatically 
insulated system in terms of fluctuations will be found (1°11) by integrating 
(11.30) over the volume of the system, using the balance equations for mass, 
momentum and entropy and applying the boundary conditions for the insul- 
ation of the system. 

(*) The use of either expression is of course correct when no relation exists bet- 
ween the time derivative of a state variable and other state variables. When this 
is not the case, however, one can have, in general, different results from the two 
expressions. In the present case it can be seen from (II.4) and (II.11) that relations 
of the mentioned kind exist. Nevertheless, since these are relations between the time 
derivative of a variable which is an even function of the particle velocities («type 

variable) and variables which are odd functions of the particle velocities (S-type va- 
riables), no possibility of confusion may arise: this results from the fact that entropy 
is an even function of the particle velocities, and is expressed by quadratic terms 
containing either «-type variables or p-type variables separately ((I.7)). 

(10) P. Mazur and S. R. DE Groor: Phys. Rev., 94, 224 (1954). 

(11) R. FrescHi, 8. R. DE Groor and P. Mazur: Physica, 20, 67 (1954). 
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Integration of (II.30) gives 


(11.42) | ee | f J,-grad AT + S Jy-grad 47, + IL: Grad Av + 
. k=1 


y 


+ Ai-c1 0AA/et — ci: (Av AB)}dV, 


were we have written, instead of the variables 7, hes A, v and i their fluctu- 
ations around their equilibrium values (*) and where (II.11) has been used. 
Partial integration of the first three terms on the right hand side of (11.42), 
with the boundary conditions of vanishing heat flow and material flows, gives 


k=1 


(11.43) | TodV= | { AT div J, + > Aji, div J, + Av-Div I! — 
: ; 


— Ai-c0A4A/0t + Av-(c Ai) B) hav à 


Using (II.4) and (11.32), and neglecting third order terms, (11.43) can be re- 
written as 


(11.44) | he | { AT (848, [t + div s,Av) + ¥ Afix(0Aox/0t + 
n k=1 
Vi 14 


+ div 0, Av) — Av-DivIl' + Ai-e 1 04A/ct — Av: (ce Ai B) +a 


or, integrating again by parts and using the boundary conditions, 


(11.45) | cav- = | {AT 0As,/ct + Y Afi, 0Aox/dt + Ai-e-1 AA] — 
BA 
v p 


— Av-(s, grad AT + > o, grad A, — c1Ai\B + Div II’) \dV È 


h=1 
Applying (11.9) and the relation 


(11.46) SEO CE > Ox grad“, = grad p + oe grad m + > 0, grad w, 
“se k=1 ‘ n 


which follows from the second law, the Gibbs-Duhem relations, and (II.31) 
G ? 


el We have Put AD TT dau AAT A, AD = v —v 
and At = t — i, (Thor Ag, Vo and i, being the equilibrium values), and have AS 


plied the equilibrium conditions, which require the constaney in time of all parameters 
the uniformity of 7, and %,, and the i eae 


vanishing of v, and iy. 


oe TS 
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equation (11.45) can be written as 


k=1 


(11.47) T, | odV=— | {AT 2As,|2 + > Aji @AQ./0t + dic 0AA/ot— 
Di F 


SE Av: (oe grad Pp dl > Ox grad wi, — ot Ai\B + Div P' )} av } 
k=1 


Applying (II.14) one has 


k=1 


(IT.48) T, | nude [iar 0As,/0t + > Ai 0A0,/8— dio 1 0AA/at + 
ki y 
+ Av: dAg/et + peAv-c1 0AA/at}dV . 
Putting 
(11.49) Ag = Ag + gee JA, 


(II.48) becomes finally 


FRE) | T's | da 


F 


= [CAT aAs, [et + À Aji, 0Ag,/8 + Ai-ot 0AAJO + Av: dAg/WAV. 
fo k=1 
For an adiabatically insulated system, as no entropy is supplied to the 
system by its surroundings (cf. (IL.28)), we have, from the integral of (11.27) 
and using Gauss’ theorem, for the change of the total entropy per unit 
time dS/dt: 


(11.61) dS/dt = | ca 


Vi 


k=1 


roa re | (AT dAs,[0t-+ Ÿ Afi, CAoufèt + Aire GAAJOt + Av-24g/2t) AV . 


V 


In order to apply fluctuation theory, the system must be energetically 
insulated, i.e. not only adiabatically, and it must be remarked that this con- 
dition for a system in which electromagnetic phenomena occur requires special 
attention. As it is shown in (%) the condition of energetical insulation does 
not change however the expression (11.51). 

In fluctuation theory we strictly need dS/dt expressed as a function of 
independent state variables. The state variables As,, dex} AA and Ag in 
(11.51) however are not all independent, because there exist n +1 relations 
(in the form of integral equations). These relations are the conservation of 
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mass of the separate components and the conservation of total energy. As 
can easily be shown, these relations do not affect our results (cf. also (**)). 

It is clear that (II.51) has the form (1.19), hence the required form for 
the change of entropy of an energetically insulated system, which is appro- 
priate. for the application of fluctuation theory (cf. (1.15), (1.16)) and mi- 
croscopic reversibility (cf. (1.23)-(I.25)). The variables As, and Ao, are of the 
a-type, whereas the components of c14A and Ag are p-type variables. (AT 
and Aj, are the corresponding x-type variables; the components of Ai and 
Av are the corresponding \-type variables). 


5. — Fluctuations. 


To these variables we apply the relations (1.15) or (1.16). This gives 


(11.52) As, (r)Q(r')AT (r') = kT, Q(r') dr — r'), 

(11.53) As, (r) Q(r') Ani, (r') = 0, (l=1, 2, ..., n) 
(11.54) As, (r) Q(r) Ai, (r) = 0, (j = 1, 2, 3) 
(11.55) As, (r) Q(r') Av, (r') = 0, | (j = 1, 2, 3) 
(11.56) Ao, (r) Q(r')AT(r') = 0, (k=1, 2,..., 0) 
(11.57) Lot) Q(r Ar) = kT, da Q(r') d(r—r'),  (k,l—1,2,..,n) 
(11.58) Ao, (r) Q(r') Ai; (r') = 0, (k=1, 2, ...,;-j = 1,2, 3) 
(11.59) Ao, (r) Q(r') Av; (r') = 0, (k=1, 2,...,; 7=1, 2, 3) 
(11.60) AA(r) Qr')AT(r') = 0, (i=1, 2, 3) 
(11.61) AA (r)Q(r)Aü(r) = 0, (i=1,2,3; 1—1,2,...,n) 
(IT.62) AA (r) Q(r') Ai; (r’) = ct Q(r') Ô(r — r'), (i, j =1, 2,3) 
(11.63) AA{r) Q(r')Av;(r') = 0, (@, FET, 23) 
(11.64) AG: PM) Q(r)AT(r’) = 0, (i= 1, 2,3) 
(11.65) AG; (r) Q(r') Ar) = 0, (i= 1, 2,85 l=, 2 oa 
(11.66) AG; (r) Q(r') Ai; (r) = 0, (i, j= 1, 2, 3) 
(11.67) AG. (r) Ar) Av;(r) = kn 6521!) de — r'), (i, j= 1, 2, 3) 
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where k, | indicate chemical components and à, j Cartesian components. Q(r) 
is an operator of the form (1.17). 

The formulae (IL.52)-(11.67) will be used in the calculation of the reci- 
procal relations between phenomenological coefficients. 


6. — Microscopie Time Reversibility. 


The equations of motion of the individual particles are invariant under 
time reversal, when also the external magnetic induction B,(r) is reversed in 
every point. This property can be expressed with correlation functions by the 
following expressions, which are examples of relations (I.23)-(1.25): 


(11.68) ss, (r) (0/01)4s, (r') (By, Bo} = As, (r') (0/0t)As, (r) {— By, — Ba} , 
(11.69) As, (r) (0/00) A0x(r'){Bi,B} = ~— Ae:(r’) (0/00) As, (r) {— By, — Bo}; 
(11.70) Zoe) 000407) Bo, Bo} = Aer) (2/t) Age (r) {— Bo, — Bi}; 
(1.71) At) (@/e)AA,(r') (Bo, By} =— 441) (0/61) As, (r) {— Bo, — Bi}, 
(11.72) Desir \(O/eAAAr’) (Bo, Bi) =— 14,7) (8/et) dae (r) Bo, — Bah, - 
(11.73) At) GAA) Ba, By = AA) (0/01) Ar) {— Bo, — Bo} » 
ari) Agr) COAG MB, By = (1) (0/4) (Bo, Bd, 
Mr) Asse) [eH AG: (1) {Bo, Bi} = — AG; ) (0/45, (r) {— Bo, — Bo} , 
(1.16) dont) (010047, (0) Bi, Boy =— AG ©) 0A) Bo, Ba, 
(1.77) Zend F') {Bo By = 4G, (FV (B/C) AAMr) (— By, — Bo. 
with k,1=1,2,...,n and i,j =1,2,3. The averages in the members of the 


equations (11.71), (11.72), (11.75) and (11.76) have opposite signs, the reason 
for this being, that As, and Ag, are even functions of the particle velocities, 
whereas JA and Ag are odd functions of these quantities (cf. chapter I, $$ 3, 4). 


7. — Regression of Fluctuations. 


In order to apply fluctuation theory, it is necessary to derive expressions 
for the decay of fluctuations, i.e. to express the time derivatives of the state 
variables As,, Jor, AA and Ag, appearing in (II.51), as functions of the 
conjugate variables AT, 4fix, Ai and Av from (11.51). 

In order to achieve this for 24), we cannot just use one of the equations 
(11.34) or (11.39), since this would not give 24/0 as a function of the conjugate 
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variables mentioned. We can, however, proceed as follows. Introducing (IT.40) 


into formula (II.11), we get 
(11.78) i= > (ex — en) Jr. 


Inserting (11.39) into (11.78) and then solving c10A/ot, we find 


n-1 
(II.79) ¢10A/ét =— A,,:grad T— YAu-grad (#:— Un) — A,,-t + o1(v AB), 
1=1 
with the abbreviations 
nes 
(11.80) Ae == (Cx: Tr Cn) (Cx FE Cn) Lira ’ 
be’ kk" =1 
(11.81) AN = (ee TS AN LL VER) 
= 
(II.82) Way Se (ee TS CAT Le: . (0) = 1, 23 very RL) 
k= 


(A. in (II.80) means the inverse tensor of A,,). So, in this way 04/0 has 
been found as a linear function of the proper variables. 
Substituting (II.79) into (11.38) and (11.39), we obtain 


n_-1 
(11.83) Je Fui À. grad 1 DA, ‘grad (lu a Ha) ole ape i ’ 
1=1 
n-1 
(11.34) J, = —A;j,:grad TS A,,-grad (Br) + Aye t, (b=1,2,, 1 0 
v1 
with the abbreviations 
n-1 
(11.85) Ie = Les x (Cx TE: En)(Exr — En) Lise EE Ly; o) 
PORCA 
n-1 
(11.86) A. = Li la DI (ex nr Cn) (Cx» xi En) Lis "AS i Lin ’ (l — 199 = 1) 
ie! ke" =1 st 
n—1 
(IL.87) NG == (ex = CALE A, 5 
k'=1 
n-1 
ESSA ee (a PI Ne 
) k ks sz Cn) (Cx» Cn) Lin A? Li 6) (k = 15 2, Ve 1) 
n-1 
cote = LY (et eee | 
) kl kl = 2 (6x En)(Exr En) Lx: Nr 3 (k, [= di D cog LES 1) 
En 


(IL.90) Ax == (ex: LE Cn) Lx: DA 


Asi i 
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a 


We need to express the tensors L explicitly in terms of the tensors A. From 


(11.81), (11.82) and (II.85)-(II.90) we easily find 


(11.91) AN AA, 

(11.92) DNA han, (ORO A) 
(11.93) De A, A ASA, (hE 1, 25.501) 
(11.94) TESINA SATIN DONA (e 


Since we have n? tensors L (i.e. 9n? coefficients) and (n+1)? tensors A (12906 
9(n +1)? coefficients), there must exist 9(2n+1) relations between the As. 

We can now write the equations for the decay of fluctuations. Inserting 
(11.83) into (11.32) and neglecting second order terms in fluctuations, (which 
is in agreement with the fact, that we have neglected third order terms in 
equation (11.51)), we obtain 


nol 
(IT.95) dAs,/et = div (A,-grad AT) + Y div {Ay-grad (4%,— Apn)} — 


1=1 


— div (A,,- Ai) — div (s,” Av), 


where s is the equilibrium value of s,. 
Introducing (11.84) into (11.4) (for k= 1, 2,-..; n— 1), we find 


(11.96) 04 or/0t = div (A,,-grad AT) + > div {Ay-grad (Aji, — 4Un} — 


1=1 


— div (A: Ai) — div (of Av), Eee E ed) 


with o the equilibrium value of gx. From (II.4), (11.40) and (II.34) we can 
also find an equation for d4g,/0t as a function of the variables AT, Afii— An, 
Ai and Av. As a calculation however shows, this equation will not provide 
any new results in the derivation of reciprocal relations (see also § 8), and 
therefore we have not written it explicitly. 

Equation (11.79), written in terms of fluctuations, becomes 


n-1 


(11.97) ec 0AA/et = — À, grad AT — > Au grad (Aji — Afin) — 


1=1 


A, Ai +e (Av) B,), 


22 R. FIESCHI | 


B, being the external magnetic induction, which is supposed to be independent 


of the time. 
Finally we find from (II.14), with (11.9), (11.41), (11.46) and (11.49) 


(11.98) 0AG,/ ot = — so oA T 0x; Es > 0), 9 OA in| OX, a cA(ATA B,); + 


k=1 


+ one (Av) B,); al 


+ D (0/00, sia L(0A0n/0%m + 0AUm|0%n 105 (p= 1535S) 
u=1 mym=1 
where 0 and e® are the total density and the total charge per unit mass 
respectively, in the state of thermodynamic equilibrium, and where à = 1, 2, 
3 denote Cartesian components. 
Equations (II.95), (11.96), (11.97) and (II.98) are the desired equations, 
which express the time derivatives of As,, 4o,, AA and Ag in terms of AT, 
Afi,, Ai and Av. 


Remark. In general the tensors L, A and the coefficients Lima Will 
depend on r and B*(r), where B*(r) is the magnetic induction at the point r, 
in a coordinate system, which moves with the same velocity as the (average) 
barycentric velocity v in that point. In fluctuation theory, however, the 
(average) barycentric velocity is zero, and therefore B*(r) becomes B,(r). 


8. — Reciprocal Relations for Heat Conduction, Diffusion and Cross-effects. 


The reciprocal relations can now easily be obtained by the preceding results 
(cf. chapter 1, § 5-6). One should insert the time derivatives of the state 
variables, expressed by the phenomenological laws (I1.95)-(I1.98), into the equat- 
ions which express the macroscopic reversibility, and then use the results of 
the fluctuation theory. 

We shall now derive reciprocal relations from the expressions (IL.68)-(11.73) 
for microscopic reversibility. Substitution of (11.95) and (II.96) into (IL.68), 
(II.69), and (11.70) gives: 


ey , ! i ! Le) ° I 1 ~ ~ 
(11.99) As,[div'{A,(B,)-grad’AT’} + > div'{A,(B}) : grad’! (Af — Ap!) — 
1=1 
— div'{A,(B;): Ai} — div'(s® Av')] = 
el Sie wer Ce te 
= As|div{A,,(—B,)-gradAT}+ Y div{A,,(—B,)-grad (A Li—AU,)} — 
t=1 


— div{A,,(— B,)- Ai}— div(sAv)] , 


' 
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A . I I fi I ! Er . Bae? i I men, ler] 
(11.100) As, [div'{A,,(B;)-grad'AT'}+ ¥ div'{A/,(B}) grad! (An) — A401)} — 
1=1 


AS div'{A,, (B;): Ai - div'(o{ Av')] = 


! . c n=1 
== Ao, [div{A,.(—B,)-grad AT} +> div{A,(--B,)-grad(4%- 47,)} — 
l=1 


— div {A,,(— B,)- Ai} — div (s® Av)], (ZERO 


CA! 1 r 1 TE 1 I 1 1 
(11.101) Ao,[div’ {A,,(B,)-grad’ AT + 2, div {A ,.(B,) : grad'(Au,, —- Au,)} — 


— div’ {A:.(B:): Ai} — div'(o{®Av')] = 


n—1 
= Ag;[div{A,(--B;)-grad AT} + > div{Auw(—B):grad (Ar —Aû,)}— 
k= 


— div {A,.(— Bj) Ai} — div (0 Av)], Get) ee 


where dashes indicate dependence on r’. 

We also could have taken (11.68), (11.69) and (11.70) for & and (or) { equal 
to n. Then we should have to substitute an expression for 0/0,/0t similar 
to (11.96). As it is already remarked in § 7, we do not get new results in this 
wav. Therefore we have only taken formulae with k, {= 1, 2,.... n —1. (For 
the same reason this is done below). 

By means of (I1.52)-(I1.59), the equations (II.99)-(II.101) can be written as 


(11.102) div’ {A (B,) -grad' d(r — r')} == div {A,,(— B,)-grad d(r—r')}, 


n-1 
(11.103) div’ {A{(B;)-grad' 6(r — r')} = div{ D Au(—- Bo) (du — db) 
t=1 
«grad d(r— r')}, (k= 1,2, Sn 1) 


(II.104) div’{ x A',(B,) ‘(a — Onn) grad’ d(r — r')} + 


1'=1 


n—1 
= div{ > Mie(— Bo): (Om — Om) grad d(r — rey ae I al 


k'=1 
The d-functions can be eliminated by multiplying with an arbitrary funct- 
ion f(r) and integrating over r. From (II.102) we obtain 


/ 


(11.105) div {At (B,)-grad f(r)} = div {A,,(— B,)-grad f(r)} , 
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where Kronecker 6’s have also been eliminated, and the left hand side of the 
equation has been twice partially integrated. The symbol + indicates trans- 
posing of the Cartesian components. Since f(r) is an arbitrary function, we 
obtain, equating the coefficients of the same first order derivative of f(r), 


(11.106) Div A,,(B,) = Div Ai, Bi) - 


Moreover, equating in (11.105) the coefficients of the same second order deri- 
vative of f(r) we have 


(11.107) AB) +A5.(B) =A,(-B) +A;.(- Bi). 
In an analogous way we obtain, from (II.103) and (II.104), 


| A,;(B,) + At (B,) À; B,) 4 Ai ( B,) and 
(11.108) | 


Div A,,(B,) = Div A‘,(— B,), (k =1,2;.2, 2— 1D} 


| Ax:(Bo) + Aj, (Bo) = A,(— B,) A AGA B,) and 
(11.109) i 
DivAy(B,)= DivAi,(— B,) . (k,t=31, 2-00 =) 


Substitution of (11.95)-(11.97) into (II.71)-(II.73) gives: 


: ‘ ; n_1 
(11.110) cAs,[—A,(B,) grad'AT'— > A, (B;):grad' (Au, — Au,)— 
1=1 ie 


— Aj,(Bi): Ai’ + e(Av' \ B))|, =— AA div {A,,(— B,)- grad AT} + 


TEEN EL rr 
ae 2 diviAn(— By) «grad(40,—42,)}— div{A,(- B,)- Ai} —div(s Av), 
(j e 1, 2, 3) 
CAES ln ¢Ao,,—A/,(B,)- grad’ AT'-- SA! (Bi) “grad (Ap Ae 2 
1=1 i ia 


— A..(B) di + e-(Av' \ B,)|,;==—AA,(r')[ div {A,..(—B,) grad AT} + 


n—1 


== 2 div {Ax (—B,) . grad (At.—An,)\— div{A,.(—B,) È A De div (0% Av), 
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Pa PATE n ie Sh ee 
GERS) cAA—A,,(B;) - grad'AT'— ¥ A/(B,)-grad! (Au, — 4) — 
L= 1 


— À,,(B;):Ai + c(Av' \ B))); = c4Aj[— A..(— Bo) - grad AT — 


— ¥A,(—B,)-grad (Aji, — AË,)— A, 


t=1 


(-B):di—-c‘A4vANBò)l, 


e 


(i, 1 = 1, 2, 3) 


where the indices à and j at square brackets denote the components of the 
vector between the brackets. 
The equations (11.110)-(11.112) can be written, with (I1.52)—(I1.63), as 


! 


(B,)(8/2a,) ô(r— r') = 


es ju 


(11-113) >A 


3 
= È (0/0x,) {A se,us(— B,) di; or r')} ’ (j == 1, 2, 3) 
i aay I I / 
(11.114) S Al, 5, BI n Sen)(8/a",) der) = 
1=1lu=1 


3 3 
(11.115) D AL B,) Oi 0(r—r') = D ANSE B,) Ojn O(r mL (id =1,2, 3) i 


where i, j and w denote Cartesian components. Eliminating again Kronecker 


ds and 6-functions, we get 

(11.116) AB ABBA Be), 
CLE?) Ara Bo) = Hit Bien Obs AB) At (—B,), 
(11.118) dBA Bio NB At (— B,), 


Witkin}; Du, 2, sand b= 1; De Cle 

The relations (11.116) and: (11.117) are examples of reciprocal relations, 
which follow from expressions for the microscopic time reversibility relations 
containing both even and odd variables (7: 

The reciprocal relations (I1.107)-(11.109) and (11.116)-(1L.118) have to be 
considered as intermediate results, because the A’s have only been introduced 
as convenient abbreviations in fluctuation theory. The desired reciprocal re- 
lations for L,,, Lu, Lis and L,, can be found from equations (11.91)-(11.94) 


spat eect 
me n 
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with (IT.107)-(11.109) and (II.116)-(11.118). From (II.91), (11.107), (11.116) 

and (II.118) follows 

(11.119) L,,(B,) + Li.(B) =[A..(By) + Al,(Bo)] +[A,.(Bo):A,,(Bo)-A..(Bo) + 
+ {A,.(B,) A; (By) -A.e(Bo)}"] = [AS Bo) + Ase(— Bo)] + 
LEA (Bi) A" (Bi) AB) +A1,( By) Az." (Bo) -Aj,( Bo) = 
=[At,(—B,) +A,,(—By)]+[ {As.(— Bo) -Az(— By) -A.s(— B)} + 


+ AE B,) AL B,) Ni B,;)] == L,.(— B,) + Li (— B,) ’ 


e 


and, in an analogous way, 
(11.120) Div L,,(B,) = Div L',(— B,) . 


In the same way it is possible to derive from (II.92), (II.93), (11.108) and 
(IT.116)-(IT.118) 


{ L,.(Bo) + L',(B;) = | om B,) + Li,(— B,) and 
TL.121 
| Div L,,(B;) = Div Li,(—B,) , (k =1, 2,..., n—1) 


and from (11.94), (11.109), (11.117) and (11.118) 


L,,.(B,) + Li,(B,) =Ly(— B,) + Li(— B, and 
11.122 | 
i | Div Lu(B,) = Div Li (— B,). (ic Lae 2 oe 


With the convention of taking all phenomenological coefficients zero in empty 
space, and the fact that these coefficients do not depend on the shape of the 
sample it follows from (IT.119)-(I1.122) 


(11.123) L,,(B,) = Lt,(— B;), 
(11.124) L.;(B) = Li.(-B), (K=1,2,..., n—1) 
(11.125) L,.( By) = Li.(— B,) . (k,l — 1,2, +» N—1) 


Taking into account (II.36) and (II.37), relations (11.124) and (II.125) are also 
valid for k and {, or both, equal to n. 

Equation (11.123)-(11.125) are the reciprocal relations for heat conduction, 
diffusion and cross-effects in anisotropie media in the presence of 


an electro- 
magnetic field. 


hh E 
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9. — Reciprocal Relations for Viscosity. 


Here we shall consider the microscopic time reversibilities (II.74)-(11.77). 
Inserting equation (II.98) into (II.74), and using (11:64)-(11.66), we obtain 


(11.126) AG;| oec-(Av' \ By); + 


Fr ee ee in tn ip ME LS Ci 
+ D (0/02) { D Lis cnn(Bo) $ (0A0n/d0m + 040m/ 0%) $] = 
m,n=1 


u=1 


RG date GIORIO 3 
= Ag, [— bee (Av AB): HE > (9/22) { > Lei mn(— B,) 
3 (04 Dr] Cm At CAVm| 2%») }] . (2, j=1, 2, 3) 


With (II.67) we find for (11.126) 


(11.127) > (d/dx,) SD ae? (0/00, ) tr — r') + bim( 0/6x,)d(r — r')} | = 


u=l mn=1 

3 3 

=> (0 OD) [ > L vai B,) {0 nl 0] 0%m)O(T—r") "i Ò;ml 0/0xn)d(r — r')}|, 
u=l mn=1 : 


(i, j — 1, 2, 3) 


since a calculation shows that 
(11.128) Ago ec (Av! B,), =— AGj0ec-(Av \ By): 


is an identity, when (II.67) is applied. Eliminating 6-functions and Kronecker 
6’s we find from (II.127) 


(11.129) L, venia) sin Lénine Bo) = Le mal B,) + Len en B,) 9 


3 


3 
(11.130) 2 0/ OX.) Las (Bi) = > (0/0Xm) Lens an B,) ? 


= m=1 


with i, j, m, vw —1,2,3. With the usual boundary conditions (see § 8), we get 
(11.131) L, us), (me) (Bo) a Lente B,) . (4, Wis mM, U = its 2, DI) 


These are 15 reciprocal relations amongst the 36 viscosity coefficients. This 
leaves us with 21 independent coefficients. 

When (II.95)-(II.98) are inserted into the equations (II.75)-(1I.77) and 
when (II.52)-(II.67) are applied, we find three identities for (11.75), (11.76) 
and (11.77). As there arise no special difficulties, we shall not give an explicit 


calculation here. 
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10. — Case of Non-Vanishing Average Barycentric Motion. 


So far we have found only reciprocal relations (equations (11.123)-(11.125) 


and (II.131)) for the case, in which the (average) barycentric motion is zero. | 
When, however, the centre of gravity has an (average) velocity different from 


1 + iq P * 
zero, these reciprocal relations still hold, when B, is replaced by B*, i.e. the 
magnetic induction in a coordinate frame moving with the centre of gravity 
(cf. remark at the end of § 7). So we find 


(11132) L,,(B*) = L' (— B*), 

(11.133) L,,(B*) = Li (— B*), (k =1, 2,.... n—1) 
(11.134) L,.(B*) =- Li, B*), (k,1=—1,2,.., n—1) 
(11.135) Laino Bi Lee DE (ij, mu=1,2,3) 


These final results are the reciprocal relations for heat conduction (11.132), 
diffusion (II.134) and cross-effects (II.133) and also for viscosity (11.135) in 
mixtures of charged components with arbitrary movements in an electro- 
magnetic field. i 


11. — Anisotropie Mixtures of Uncharged Components. 


From the general case which has been treated in the foregoing one can 
easily derive the more particular cases in which the physical situation implies « 


that not all the fluxes and forces are present. 

The case of the anisotropic mixtures of n uncharged components, in which 
the phenomena of heat conduction, diffusion and cross-effects, and of viscous 
flow may occur, has been explicitly treated by DE Groot and Mazur (5) as 
a first example of their generalized Onsager theory. If we simply put e, = 0, 
the phenomenological equations (IL.38), (11.39) become 


n—1 
(11.136) J, = — L,. grad T — > L,,: grad (a, — Hu) 
1=1 
m— 1. 
(11.137) Je — L,.: grad = D: L,,-crad (Hi an Un) , (RES DI tees n—1) 
t=1 
where now #,=w, + w, (the system may be subject to external conservative 


forces). The equations for the viscous flow already coincide with the corres- 
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ponding equations which those authors have given. The reciprocal relations 


for the phenomenological coefficients, which have already been demonstrated 
in the previous sections, still hold. 


12. — Galvanomagnetic and Thermomagnetic Phenomena in Anisotropic Metals. 


Let us consider a system, consisting of a rigid ion lattice and of electrons, 
in an electromagnetic field and with a non uniform temperature. Equations 
(11.38), (11.39) reduce in this case, taking the barycentrie velocity zero, to 


(11.138) J= = Ly grad T— L,,' {grad (a, = i.) — (e, — ei)e ? 0A/ct}, 


é;)c-! 0A/0t}; 


(11.139) J,= — L,,-grad T= L,,: {grad(f, — fi) — (e. 
moreover (11.36) reads 
(11.140) J,=—J.. 


(index e indicating electrons, index è ions). By introducing the electric 
current and with (11.140) we obtain 


(II.141) i + CeOeVe + C:0;Vi = (e, “sd eid o 


Using now the approximation ¢,>e, (the charge density per unit mass of electrons 
is much larger than the charge density per unit mass of the ions), we can 


finally write 


(IT.142) JL" vrad TTL. "{orad (4 Hi) 40° 0A Ct} . 


(11.143) eo ordi — 6," (grad (ue pi) — 6 * 0A / et} . 


Reciprocal relations (I1.123)-(11.125) reduce to 
(11.144) L;.(B)=LL(— Bo), L..(B;}= Li Boysen La(B)=L Bo): 


These equations enable us to treat the galvanomagnetie and thermomagnetie 
phenomena in anisotropic metals. They will be reduced in the next chapter 
to the case of isotropic metals, in which situation galvano- and thermomagnetic 


phenomena are actually studied. 


3 Supplemento al Nuovo Cimento. 
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Remark. These phenomena have been already treated, from the point of 
view of thermodynamics of irreversible processes, by CALLEN (2%) and by 
Mazur and PRIGOGINE (4). Those authors, however, do not give a proof of 
the reciprocal relations: as it was usually done, they simply suppose that 
those relations are examples of the relations which ONSAGER (*) derived for 
scalar processes from the property of microscopic reversibility. 

The same case of an anisotropic metal in a magnetic field has been expli- 
citly studied (1') by introducing the physical simplifications from the beginning 
in the fundamental laws. We should like to say, however, that the reduction 
of the general case to a more particular case, and the introduction of the 
simplifications in the phenomenological equations, allow sometimes a more 
complete description of the phenomena. 

The equations of this section, as it may be easily seen, include the case 
of electric conduction in anisotropic crystals in a magnetic field at uniform 
temperature, which has been explicitly treated by Mazur and DE GROOT (1!) 
by using their extension of Onsager’s theory of reciprocal relations. 


CHAPTER IIT. 


EXPLICIT EXPRESSIONS FOR THE MEASURABLE EFFECTS 
IN ISOTROPIC METALS 


1. — Introduction. 


In the preceding chapter we have derived the Onsager relations for heat 
conduction, electrical conduction and cross-effects in anisotropic metals in the 
presence of a magnetic field. The purpose of this chapter is to derive the 
expressions for the observable coefficients usually defined tor the case of iso- 
tropic metals (15). 

The coefficients for the galvanomagnetic and thermomagnetic effects have 
already been given in local form (!?41), i.e. the various relations between 
electric field, temperature gradient, heat and electrical flow are expressed in 
terms of the phenomenological coefficients at a certain point of the system. 


PSP Carn? Phys: ev. 73 01310 (1948). 
ON Tal. ans Phys. Rev., 85, 16 (1952). 
Ps 
R. 


Mazur and I. PRIGOGINE: Journ. Phys. et Radium, 12, 616 (LOST): 


FiEscHI, 8. R. DE GROOT and P. Mazur: Physica, 20, 259 (1954). 


15 


ds) 
VEN) 
ta 
i) 
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However these local coefficients do not necessarily correspond to the experi- 
mentally measured coefficients (1518), because « heterogeneous and homogeneous 
thermoelectric effects » and Peltier effects may also play a role, in addition to 
the purely galvano- and thermomagnetic effects. Consequently, of course, the 
relations among the various effects have also to be modified. 

Besides the effects which have already been studied in local form by other 
authors (1215), mentioned above, we have considered also the analogous effects 
in a longitudinal magnetic field, and the Peltier effect in a transversal and 
longitudinal magnetic field (21819), In particular we remark that Thomson’s 
second relation is still valid in the presence of a magnetic field. 


2. — Phenomenological Equations. 


Let us consider an isotropic metal in a magnetic field B, parallel to the 
z-axis. As can easily be seen, the phenomenological equations (11.142), (11.143) 
for this case, neglecting the time derivative of the vector potential 04/0, 
(i.e. assuming that the electron motion is sufficiently slow to allow us to neglect 
the effect of an induced magnetic field), can be written as (*) 


WD) bp = — Lee: O(e— Bs) OB — Leo aye U_U) — Lera 0T/00— Les ny OT | OY, 
+2) Vy == LD AU Hi) [Oa — L ee, ra? e 10( Me a U:)/0y =F dolce Lice zeoL | CY; 
Bt, = — Lit (U_U) Losa OT [02 , 
I ge ex O(U— [;)/0a — Le ryes U_U) OY — Lise rnO Lf 00 — Liss my OL} CY, 
.D) di aa Lia o(U— n.) | 0a — Lens 30 fA) | OY ao Liss OL Ce Leo Do 
BI, = — Lic Bi) 0% — Lu oT | dz , 

(6) W. MEISSNER: Handbuch der Experimental-Physik, 11 (Leipzig, 1935). 

(7) W. GerLacH: Handbuch der Physik, 13, Chapt. 6 (Berlin, 1928). 

‘ (18) L. L. CAampBELL: Galvanomagnetic and Lhegmomagnete effects (London, 1923). 
ee G. LASKI: Hanabuch: der mo 13, re Alina no 


of equations ai. 1- an; 
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The isotropy of the system requires, in addition to the following relations 


between the coefficients: 


| De LE Le YU 3 Lo, ax = Lig uy 9 Lie, ax — Le, ? Lg xx — Lis ’ 


LT 
I L L =—L L —L 
vu — LO yxy +/es,ny — — Mes, yuo ““se,ay 80,0) ~~ 88, bY TT 88, V2) 


ee, LY 


which have already been taken into account in the equations (IIT.1)-(III.6), 
also that the coefficients L.. 22; Los say Des arr Lineo) Lee ar Vea ves Liscia, ond! Dio 
are even functions of the magnetic field, and that Les, Desovy Ls, and 
Lm are odd functions of the magnetic field. 

Moreover the Onsager relations (11.144) hold between the coefficients: 


(8). ZL. o(B)= La, B)p L..,4(B) = Lau B), La B= Be 


where i and j denote Cartesian components.’ Some of these relations already 
follow from the isotropy of the system, the really new relations being 


Vp add sy) = re B) ’ 
or, since Ls; .. is an even function of B, 


(IIL.9) DEB 


Joy; oy(B) = boe B) ’ 
or, since Le» is an odd function of B, and from (III. 7), 


(11.10) Lie ey( B) = Lo, on B); 


Tee BD) = Ly al— B) 9 
or, since L,,,. is an even function of B, 


(111.10a) L,s,:(B) = Le dey) 0 


Remark. Other authors (12141) have studied the case with fluxes and forces 
limited to the «-y coordinate plane. We have considered the more general 
case of fluxes and forces of arbitrary direction, in order to be able to deal 
also with the effects in a longitudinal magnetic field. The isotropy of the 
system requires that the coefficients with second sets of indices 20, ZY, xe 
and yz vanish. So the equations (III.1), (III.2), (111.4) and (IIL.5) coincide 
with the equations which the quoted authors wrote down. In addition to these 
we have the equations (III.3) and (III.6), which are independent of the fluxes 
and forces in the x-y plane: they are the same equations which describe 
the thermoelectric effects, but the phenomenological coefficients are here of 
course functions of the magnetic field. These equations will enable us to treat 
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the effects in a longitudinal magnetic field. (In the anisotropic case the effects 
in a longitudinal field contain of course also cross-terms, pertaining to the 
perpendicular directions). 

Mazur and PRIGOGINE (!*) pointed out that a different form of the pheno- 
menological equations provides a more convenient way to deal with the 
problem in the case of a metal. Their form differs in two respects from the 
form used above and in (12): 1st by choosing as independent variables the 
quantities which are directly controlled in the experiments (electric currents 
and temperature gradients), i.e. by changing the role of the electric flux and 
force, it is possible to obtain very simply the explicit expressions of the effects 
and the relations among them as functions of the phenomenological coefficients ; 
2nd moreover, a different choice of the fluxes and forces explained hereafter 
allows also the treatment of the effects in which the electric current does not 
vanish in the direction of the «adiabatical » insulation (5): 

The new set of variables is obtained from the former set by à linear trans- 
formation which leaves the entropy production o invariant. The entropy 
production (11.37) in our case is | 


3 
(111.11) o=—T1Y [J .< el | ex,)— J ¢j(O(fe— fi) | 0a;) ]. 
j=1 
The phenomenological equations (III.1)-(IIT.6), where è; and J,; are the flows, 
— e! a,— R)%%; and 0T/èx; the corresponding forces, are based on this 
expression. With the linear transformations 


(III.12) LETI), 


where J; is called « reduced heat flow », s, and s; are the specific entropies 
of electrons and ions, (III.11) can be written as 


3 
MS or JE) Ie e;) d|dx; — 
j=1 
cet ou ps3) | 020; jp conet}) , 


or, since the total pressure and the concentrations of electrons and ions can 
be supposed uniform, ((grad /4)7,p_ const — 0, (grad p:)ppoconst = 0) and with 
(11.141) and the approximation @, > @;, 


3 
(III.14) o=T-+Y{-J_,T-(0T/dw;) + i;E;}, 
j=1 
where the electric field E =— grad has been introduced. 
It is now possible to write linear phenomenological relations which express 


(*) In non-metallic systems the interchange of the role of electric flux and force 
is also preferable; but the different set of fluxes and forces is not suitable (15). 
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the fluxes i; and J-, as functions of the forces — T-10T/0x; and H;. As it 
has been said, however, it is convenient to express the variables J;; and E; 
in terms of the other quantities. 


IFPI) Joe Lig Lug La la eee 
Hide) Bon a, CL TIA 
oO IERI 
(LILAS) Fag Lge ect bt Lu oi Dsl dle La D OT [Ou 
DI SJ Lago Less + Laney OT] Cn Else CRAQUE 
(111-20) © Sg = Dt, 1D, 7 T eT yor 


Here the isotropy of the system has already been taken into account; from 
this .condition, it follows: also tbatL, nel, Lita a sages LR 
and L,,., are even functions of the magnetic field, while De. ay, Lesevy Lse,av 
and L,;.., are odd functions of the magnetic field. The validity of the Onsager 
relations is conserved by a linear transformation of fluxes and forces (1). Ex- 
cluding the relations which are already given by the isotropy of the system, 
they read 


(111.21) LiBeR 
(IIL.22) LAB RL ee BEL B); 
(111.23) Dre (B) a ep iy Lee (= B) a Fa LE Zz (B)} ’ 


where the minus sign results from the interchange between the electrical flux 
and force. 

With the aid of these relations the independent phenomenological coef- 
ficients reduce to six for the equations of the components in the «-y plane, 
perpendicular to the magnetic field: Lx», Le an Les Dal ANIELLO 
and to three for the equations of the components in the direction of the mag- 
méticielden 21 o ands Ll 


88,22 * 


3. — Local Effects and Relations among Them. 


The effects we study are defined in MEISSNER (15) (see also GERLACH (1°)). 
The galvanomagnetic effects are caused by an electric current; the thermo- 
magnetic effects are caused by a heat flow. In the transversal effects the 
primary current (heat or electric current) is perpendicular to the produced 
effect; in the longitudinal effects both variables have the same direction. An 
isothermal effect is defined in such a way that the temperature gradient in 
the direction perpendicular to the primary current is zero; for an adiabatie 
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effect no heat is allowed to flow in the direction perpendicular to the primary 
current. 

Symmetry arguments on the reversal of the magnetic field suggest that 
the transversal effects are in first approximation proportional to B, the lon- 
gitudinal effects to B®. Several authors (#77) therefore divide the coefficients 
of the various effects (see following table) by B or B? respectively, in order 
to obtain new coefficients, which are «independent » of B. We prefer, how- 
ever, not to express explicitly the dependence of the coefficients on the mag- 
netic field, since the experimental deviations from the foregoing simple de- 
pendence are too relevant. The even or odd parity of the effects with respect 
to the magnetic field is already contained in the even or odd parity of the 
phenomenological coefficients, as is shown in the preceding section. 

Below we give a table with the definitions of the various local effects (14, 
with a transversal magnetic field, of the theoretically possible 8!/2!3!31=—560, 
and 4 with a longitudinal magnetic field), and their local expressions in terms 
of the phenomenological coefficients (12:14), which are easily derived from the 
equations (III.15)-(1II.20). 

When the isotropy of the system is taken into account, the 14 effects with 
the transversal field are expressed in terms of 9 independent phenomenological 
coefficients (cf. table); hence 5 relations exist, connecting the effects with each 
other, which are simply due to the symmetry of the system: 


(III.24) | R—R=—QP', 
(111.25) À, — À = (SE 
(111.26) (ee i =e aoe 
| | (Heulinger relations) 
(111.27) Qi —Q; = IM A 
(111.28) dg Br 


Moreover, the Onsager relations (III.21) and (111.22) reduce the number of the 
independent coefficients to 6; we then have among the effects 3 more relations: 


(111.29) A;P'= TQ°, (Bridgman relation) 
(IIT.30) api IO 
(111.31) APL HORS 


Owing to the Onsager relation (III.23), the 4 effects with the longitudinal 


field are expressed in terms of 3 independent phenomenological coefficients ; 


hence one relation exists among them: 


(111.32) APE — TQ". 
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4, — Measurable Effects. 


Let us consider a system consisting of a metallic sample A and of the me- 
tallic wires B, which connect two corresponding points on opposite faces of 
the sample with a condenser © (fig. 1). We call: (Ag); = (pi 92); the potential 

difference which is measured (for 
metal 8 instance) with a condenser connected 
with the corresponding points on the 
faces parallel to the plane «,;2, (4, =a, 
Wy = Yt, = 23.1,), b= 1,2, S347 
(the figure illustrates the case 1=y); 
(AT), the temperature difference between 
ne BB, SFE the corresponding points on the faces 
3 parallel to the plane w,x7,(i4j <k); I; the 
total electric current which is measur- 
ed in the direction #,; J;; the total « re- 
duced heat flow » which flows in the direction «;, Moreover we suppose that ?,, 
J-;, 0p/0æ; and 0T/0x%;, are uniform inside the metal A. Then the following 
relations hold between the above defined measurable quantities and the cor- 
responding «local» variables: 


Fig. 1. Sample with connecting wires. 


(111.33) (qu —@ 2): = | (2/Om,) da, = Aa, (dp/2a,) , 
re 
(III.34) (AT); = | (21/0) da, = Ax, (0T/ex,) , 
Az; 
(111.35) I,=| ficdo,dm= Ao, Ao’, (Aj #4) 
Ax; Ary, 
(111.36) 1, =] |J,.dv;,dm= Av, An... (i Aj Fk) 
At; Axy, 


In the experimental apparatus one has, besides the sample A, which is 
placed in the magnetic field B, wires of a different metal B and junctions 
of the two different metals (Fig. 1: points 3 and 4). Consequently one should 
pay attention to the thermoelectric effects which may be superimposed on 
the pure galvano- and thermomagnetic effects. 

When the temperature in the wires and at the junctions is different, a con- 
denser records, in addition to the potential difference of galvano- or thermo- 
magnetic origin, a thermoelectric potential difference which is given by the 
following two effects (1). 


1) The «homogeneous thermoelectric effect », i.e. the potential difference 


THERMODYNAMICAL THEORY OF GALVANO- AND THERMOMAGNETIC PHENOMENA 39 


inside the single wires, is given by 
or RS eg 38 * 
(111.37) Pa — Py = — V5"(T.— Ti)/eT , 


where the indices a and d indicate two points of a wire, Q** is the «reduced 
heat of transfer » for the electrons in the metal B (i.e. the reduced heat trans- 
ported by the unit of mass of the electrons in the state of uniform tempe- 
rature, which can be expressed by the coefficients of the phenomenological 
equations for metal B) (*). 


2) The «heterogeneous thermoelectric effect », i.e. the potential diffe- 
rence at the junction of two different metals A and B at the temperature Le 
is given by 


(IIL.38) (Da — Pr = € (Ua — Ug)r » 


where 4, and uz are the chemical potentials of electrons in the two metals (*). 

The general equation for the resulting potential is obtained by adding up 
the homogeneous and heterogeneous potential differences in the circuit. One 
then obtains, for (III.33), (III.37) and (III.38), (Fig. 1), 


(111.39) (Ap): = (gr Pa): = (Yr — Pada + (Pa— Pa)r, + (Pa— Pali + 
+ (Pa — Par, + (Ps Pes OST, T.)+e tu pa)Ti+ 
à (0p/0æ,) dae; + eus Ma), + li T3) = 


Az; 


= ¢1T-19**(T,— T, + T,—T 5) +e [Un Me) OL UT Po) + 


+ (09/85) Ace = tg (Ta Ts + TT) + 
AE es — SMI, Ps) (0w/0æ;) Aa; . 


The quantities s, and s, are the specific entropies of electrons in the metals 
A and B, and in the last step of the derivation Gibbs’ equation has been used. 


(*) In the literature two sets of «heats of transfer » are used (?°): one of the sets 
consists of quantities; each of these quantities may be interpreted as the heat trans- 
ported per unit flow of matter of a certain chemical component. The other sets con- 
sists of n 1 quantities, where each quantity may be interpreted as the heat trans- 
ported per unit flow of matter of a certain chemical component relative to the bary- 
centric motion. Since in our case c, < Ci (Ce is the concentration of electrons, e, the 
concentration of ions), it is easily seen (2) that the «heats of transfer » of the latter 
set coincide with the first n — 1 heats of transfer of the first set mentioned (in our 
case n=2). Moreover, since relative and absolute flows coincide, the corresponding 
two kinds of «reduced heat of transfer » also coincide, and no further specification 
is necessary. 

(+) Here and in the following we let e, = €. 

(20) G. A. KLUITENBERG, S. R. pe Groor and P. Mazur: Physica, 19, 1079 (1953), $ 6. 
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If the temperature gradient vanishes in the direction in which the po- 
tential difference is measured, (7, = T, = T; = T;), the « homogeneous ther- 
moelectric effect » also vanishes, the « heterogeneous effects » at the two junct- 
ions are equal and with opposite sign, and one has simply 


(111.40) (Ap); = (0p/0æ,) Ax, = — E, Ax,. 


The isothermal Hall (I) and Nernst (V) effects, and the isothermal (VII) and 
adiabatic (VIII) electrical resistivity belong to this case. 

For the adiabatic Hall (II) and Nernst (VI) effects, one has 737 7,, but 
T,=T, and T,=T, (from the adiabatical insulation condition for the y 
direction). The homogeneous thermoelectric effect then vanishes and (111.39) 
becomes 


(111.41) Mo) = (ss) Eder 


For the isothermal (XIII) and adiabatic (XIV) Ettingshausen-Nernst effects 
-the homogeneous thermoelectric effect can also be present, since the adiaba- 
tical insulation in the direction in which the electric field is measured is no 
longer necessary; (ITT.41) is then again valid only for the case that the tempe- 
rature gradient is zero in the wires which connect the metal A with the con- 
denser. In the following we shall consider this case in order to have con- 
sistency in the definitions and more simplicity in the relations among the 
various effects. In § 6, however, to treat 
J AI the Peltier effect, it is convenient to con- 
a ze sider the more general case in which a 
ig GK temperature gradient in the wires is also 
present. 
If sample and wires are made of the 
Fig. 2. - Sample in magnetic field. 8®™€ metal, the heterogeneous thermo- 


ae 


(Point 3 is at the left-hand dotted electric effect is the potential difference 
line in the metal; point 4 is at the at the «junctions» between the metal 
right-hand dotted line in the metal). inside the magnetic field and the metal 

outside: s, is then the specific entropy 
of electrons for the metal inside the magnetic field, while s, is the specific 
entropy of electrons in the metal without magnetic field (see the end of § 6); 
T; and 7, are the temperature at the points where the magnetic field drops 
to zero. 

If the boundary of the magnetic field does not coincide with the junctions 
between the two metals, but part of metal À is left outside the field, one has 
4 discontinuity points for the electric potential (fig. 2: points 1, 2, 3 and 4), 
and homogeneous thermoelectric effects both in the wires B and in the parts 
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of A outside the field (fig. 2: the dotted lines indicate the boundaries of the 
magnetic field). 

It is now easy to give the relations between measurable and local effects. 
The coefficients with an upper bar indicate measurable coefficients, while 
without the bar they indicate local coefficients. 

The isothermal Hall coefficient (I) is defined by Ri == (Aq), 4z/I,, when 
1, —0, (AT),=0 and (AT),=0. From (III.35) and (III.40) one has: 
Rt = B, Ay Az/i, Ay de = R;. In a similar way it is immediately seen that 
the following observable effects coincide with the local effects. (For the sake 
of brevity we shall not write the conditions which define every observable 
effect, as they coincide with the conditions for the corresponding local effect). 
Isothermal Nernst coefficient (V): Qf =— (49), Ax/(4 T\ Aye 0 (Agee 
matter of fact temperature and potential differences in the x (longitudinal) 
direction are not measured at the ends of the sample A. In these cases, hence 
Ax is the distance between the points where the measurements are performed 
(fig. 2, points a and b)). Ettingshausen coefficient (III): P'= (AT),42/I,= P'. 
Righi-Leduc coefficient (IV): St — (AT,) Ax|(AT),Ay = S'. Isothermal (VIT) 
and adiabatic (VIII) electric resistivity: Rt =— (Ap) y Ay Az/I, Ax; R; = Bi, 
R'—R}, Isothermal (XI) and adiabatic (XII) heat conductivity: 
Mi = — I,, Ax/(AT). Ay Az; Ai dij Aa = da 
For the adiabatic Hall effect (II), as we have said in the foregoing, the 
influence of the potential difference at one junction is no longer compensated 
by the potential difference at the other junction, since the temperature at 
the two junctions is not the same. So we have, from (III.41), 


Ri =— (Ag), 42/1. = (Az/i, Ay ARTE, Ay + e*(83— 84){— (AT) jt] = 

= Ri + e*(84— 82){(AT)|tx AY}; 
or, for the definition of effect (III) in the table, 
(III.42) Ri = R° +e (si —582)P". 
In an analogous way one finds, for the adiabatic Nernst coefficient (VI), 
(111.43) di =— (Aq), def(AT). Ay = Qù + 8a — 8,)S¢, 


where the definition of effect (IV) in the table has been used. This correction 
for the measurable adiabatic Nernst coefficient has been given already by 


JAMPBELL (1°). 
For the longitudinal isothermal Ettingshausen-Nernst coefficient (XIII) 


one has 


(TII.44) Ot =— (Ap) (AT }a = Qi + 184 Sz) ; 


eve 
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and for the longitudinal adiabatic Ettingshausen-Nernst coefficient (XIV) 
(111.45) Qi = Qi + e Usa — 85) - 

As we have remarked, these last two results are valid for at. case in which 


the temperature gradient is zero in the wires. If the temperature gradient 
is non-vanishing and uniform (and so also 7, = T,), one has from (III.31) 


(111.46) Qi = Q + {e(s.—82)— eT 193}, 
and 
(111.47) Q! = Q! + (es, — 83) — eT 4Q3*} . 


The longitudinal isothermal (IX) and adiabatic (X) Nernst effects are the 
only two effects in which the electric current in the direction, in which one 
measures the temperature difference, does not vanish. One has therefore an 
evolution of heat at the junctions, which modifies the temperature difference 
due to the pure galvanomagnetic effect. However, the effect of heating at 
the junction should be corrected with suitable experimental arrangement, as 
it would be complicated to consider explicitly its influence on the expressions 
of the two galvanomagnetic coefficients. Apart from that, the experimental 
coefficients coincide with the local ones. The conditions supposed in this case 
are hardly to be obtained practically, as we must assume the possibility of 
the existence of an electric current without an accompanying «reduced heat 
flow». Nevertheless, the case is of theoretical interest. Longitudinal iso- 
thermal and adiabatic Nernst coefficients: 


PAT) AY AIC AG. EP PP RD Po 


For the longitudinal effects in a longitudinal magnetic field again the same 
considerations hold as for the corresponding longitudinal effects in a trans- 
versal field. Two of them (XV) and (XVII) are unchanged: R’# = R™ and 
48 — }'8, The thermoelectric power (XVIII) Q4=— (Ag),/(AT), is given by 


(111.48) Q= Q'+ es, — 83), 

if the temperature gradient is zero in the wires, and by 

(111.49) QE QUE {e(s,— 83) — e7T1Q**} , 

if the temperature gradient is non-vanishing and uniform. The coefficient 


P™ (XVI) is again affected by the influence of an evolution of heat at the 
junctions, which cannot be expressed explicitly. 
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5. — Relations Among the Measurable Coefficients. 


We have seen that the local expressions of the four effects II, VI, XIII 
and XIV, where the temperature gradient does not vanish in the direction in 
which the electric field is measured, must be modified, because for these cases 
a thermoelectric potential difference is superimposed on the galvano- or thermo- 
magnetic potential difference. Consequently, some of the relations among 
the various effects are also modified. 

The relations (III.24), (111.27) and (III.28) remain unmodified because 
the local and the measurable coefficients, which appear in them, coincide. 

It is easy to see that the relations (III.29), (I1I.30) and (III.31), expressed 
in the measurable coefficients, remain unchanged, although the expressions 
‘for the measurable coefficients are different from the expressions for the local 
coefficients. 

Relations (III.25) and (III.26) are no longer valid when the measurable 
coefficients are substituted for the local coefficients. The relations then become 


(111.50) A,B? = TQ, + e*(8,— 82) 
and 
| (ITI.51) A Pi = TQ) + e (sa — 823), 


or, by eliminating the specific entropies, 
(III.52) PAPE TQ 0) 


Also the relation (II1.32) for the longitudinal effects in a longitudinal field 
is not valid any more for the measurable effects. Instead of it we have 


(111.53) AUD = TQ + es, 83). 


In this way finally all relations among measurable effects have been found. 


6. — Peltier Effect in a Magnetic Field and Thomson’s second Relation. 


In this section we wish to extend the theory of thermoelectricity (1) to the 
case in which a magnetic field is present. 

In order to study the Peltier effect in a magnetic field it is convenient to 
extend the definition of «reduced heat of transfer ». In a transversal magnetic 
field we define an «isothermal reduced heat of transfer » (Q**)?= e ia): s 


7 
44 R. FIESCHI | 


with the conditions à, = 0, OT/0æ = 0, OT/0y = 0; and an «adiabatical re- 
duced heat of transfer», (Q**)!?= e“(J;,/t2)°, With the conditions è, = 0, 
AMATO — Onn trom) the equations (III.18) and (III.19) we have 
then (Q**) = e"Dse,m, 0r, With the Onsager relation (III.21), 


(111.54) (OM ehi 


(Q**)2 = EL + Les ou Lise, al Use, ae) OF, With the Onsager relations (III.21) 
and (III. 22), 


(111.55) (hea = ONES Le HO TR LE: reg bee ies un) © 


In a longitudinal magnetic field we define the reduced heat of transfer 
(Q**)8 = e-1(J [î,)!?, with the condition oT /ozg = 0. From equation (III.20) 
we have 

(Q**)'8== eDL,,..,, or, with the Onsager relation (III.23), 


(ITT.56) (CSA = Pl . 


(The coefficients here defined are, of course, functions of the magnetic field). 
One easily sees that the following local relations hold: 


KETO TN) Cet OR HO NET (QT) ea (OI 


These relations are, in a certain way, the «local» analoga of the «measur- 
able » second relation of Thomson. 

The Peltier coefficient is defined as the absorbed heat per unit flow of an 
electrical current across an isothermal junction from metal A to metal B. 
We can define two coefficients for the case of a transversal magnetic field: 
an isothermal Peltier coefficient: (77,3), with the conditions j= 0,.0T icv — 
oT/cy =0; an adiabatic Peltier coefficient: (/7,,), with the conditions i,=0, 
oT /ca = 0, J;, = 0; and a coefficient for the case of a longitudinal magnetic 
field: (Z/,,)', with the condition 0T/0z¢ = 0. | 

The evolution of heat at the junction between a metal A and a metal B 
can be found (1) as the difference between the entropy flows before and after 
the junction, multiplied by the temperature 7, since there is no entropy source 
at the surface separating the two metals. For the isothermal and adiabatic 
case in metal A inside the transversal magnetic field we have, from (IIT.12), 
(117.18) and (III.19), 


(111.58) TJ 2): a (Jia): A ETS = (ER CLs PE é H(QS")} +Tshe 


(11:59) T(J,, na == (Ja a)a Le Ts ie i (SERE ste Lit Eta ah T's) = 


ey ok) 1B + 
eo) dr Lee 
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From (III.12) and (III.20), for metal A inside a longitudinal magnetic field, 
we have 


O111.60) TJ: J, + OP syt, = (Leese + IT 8) = e-f{(Qi*)4+ Ts,}i,. 
In the wire B we have (1) 
(111.61) DIS en (0 PTE ue (j = 1, 2, 3) 


The difference between (III.61) (i = 1) and (III.58), divided by the electric 
current 7, gives 
(111.62) (LEO NOT) + Tiss sa). 


A z 


The difference between (III.61) (i — 1) and (III.59), divided by the electric 
current 7, gives 


(III.63) Ul) = “0 = (074 + T(83 — 84)} - 


a 


The difference between (III.61) (i = 3) and (III.60), divided by the electric 
| current 7, gives 


(111.64) (Ups)? =e QU — (077)? Ls a. 

From (111.46), (111.57) and (III.62) one has 

| (IIL.65) (uni =— TQ} ; 

from (III.47), (111.57) and (III.63) one has 

(11.66) ESA 

from (III.49), (III.57) and (III.64) one has 

(111.67) (I4,)2 =— TQ". 

The coefficients Q!, Q! and Q' can properly be called « thermoelectric power 
in a (transversal or longitudinal) magnetic field »; hence it is correct to say 
that the relations (III.65), (11.66) and (III.67) are extensions of Thomson’s 


second relation, to the system consisting of a metal subject to a magnetic 
field and a metal outside the magnetic field. 


4 — Supplemento al Nuovo Cimento. 
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Remark. We should like to stress again the fact that, when the whole 
system is made of the same metal, the effects between the metal inside and 
outside the magnetic field can be studied. In such a case 8, indicates the 
specific entropy of electrons for the metal inside the field, (when the metal 
is isotropic it does not matter whether the field is transversal or longitudinal) 
and s, indicates the specific entropy of electrons for the metal outside the 
magnetic field. As a matter of fact Onsager relations have been demonstrated, 


be mala 


in chapter II, for a metal without magnetization or polarization. However, # 


they can easily be extended to the present case, in which the influence of the 
induced field is negligible. As MAZUR and PRIGOGINE (8) have shown, the 
local entropy production does not change, formally, when magnetization is 
also taken into account: the chemical potential x only must be corrected with 
an additive term depending on the magnetic susceptibility of the metal and 
.on the magnetic field. Since in our case both quantities are constant, the 
Onsager relations can be derived in the way outlined at the end of (1). 


Appendix on Tensor Notation. 


We have used the Milne system of tensor notations (2?!) in this paper. The 
exterior product of an ordered pair of vectors a, b is a tensor 


Lab components: i ade. 
The divergence of a tensor and the gradient of a vector are written as 


Div T, components: (Div T);= >, (0/0%,)T x2, 


Grad a, components: (Grad a),, — (0/0x,;)a, . 
The interior product between a tensor and a vector is denoted by 


T-a, components: (T-a);= >, Ta, 


a:T, components: (a-T), == > Glee 
Finally 


T AU = > Le (DE 


is the interior product of two tensors T and U contracted twice. 


(21) L. RosENFELD: Theory of Electrons (Amsterdam. 1951). 
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Introduction. 


In various previous papers ("8), several new interpretations of the non- 
relativistic quantum theory without spin have been proposed. Although these 
new interpretations differ in various details, they all have in common that 
they explain the quantum theory in terms of continuous and causally deter- 
mined motions of various kinds of entities, such as fields and bodies, which 
are assumed to exist objectively at the microscopic level. Thus far, those 
interpretations which have been carried far enough to demonstrate in full 


(*) Now at Centro Brasileiro de Pesquisas Fisica, Avenida Wenceslau Bras, 71, 
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detail their ability to explain causally all of the features of the Schrédinger 
equation without spin have followed one of two general lines of approach. 

The first of these general lines, initiated originally by DE BROGLIE (°), and 
later carried to its logical conclusions by one of the authors of the present 
paper (*), involves the notion that the Schrédinger wave function, V(x, t), 
represents an objectively real but qualitatively new kind of field of force that 
influences the motion of a body, which latter has a well defined location, §(¢), 
varying continuously and in a causally determined way with the passage of 
time. This new field of force produces no important effects at the macroscopic 
level, but at the atomic level it is, as has been shown in various papers referred 
to above, able to explain the characteristic new quantum-mechanical properties 
of matter, which manifest themselves strongly only at this level. 

The second general type of causal explanation of the quantum theory is 
along the lines of the hydrodynamic model, proposed originally by MADE- 
“LUNG (2), and later extended by TAKABAYAST (4) and by SCHONBERG (°). In 
this model, |¥|? represents the density of a fluid, while VS/m represents its 
local stream velocity (where ‘= exp [iS/h]). In another paper (8), this 
model has been completed with the aid of the postulate that there is a stable 
particle-like inhomogeneity in the fluid that moves with the local stream velo- 
city. This inhomogeneity plays a role analogous to that played by the body 
in the interpretations initiated by DE BROGLIE, while the fluid plays a role 
analogous to that played by the Ÿ field. 
In the present paper, we shall find it convenient to work in terms of the 
hydrodynamic model, which we shall extend with the aid of various new pos- 
| tulates to the Pauli equation describing an electron with spin. We shall there- 
fore now summarize a few important aspects of the hydrodynamic model 
without spin, in order to facilitate the description of the nex features that are 
needed for the treatment of spin. We first write Schrédinger’s equation in 
the well-known form 


spec VS 

(1) om + aiv( Re |=0, 

as (VS)? ie VR 
(2) ot tam rai da 


In terms of the Madelung interpretation, equation (1) then represents the 
conservation of fluid, while (2) represents the equation which determines the 
velocity potential, S, in terms of the classical potential, V, and the « quantum 
potential» U=— (42/2m)V?R/R = (h?/4m)((Vo)?/2e? — V2o/0). As shown by 
TAKABAYASI (4) and by SCHONBERG (7), the quantum potential can be inter- 
preted in terms of a kind of internal stress in the fluid, which depends not, 
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however, on o itself as is the case with pressures and other internal stresses | 
in an ordinary macroscopic fluid, but rather, on derivatives of o. 

In one of the papers referred to previously (5), the additional assumption 
was made that the fluid is in a state of irregular fluctuation resembling turbulent 
motion, in which the actual density 0, and the actual velocity v, fluctuate 
more or less at random around |¥|? and VS/m respectively as means. In 
these fluctuations, it is permissible to assume that v has vortex motions, which 
however average out to zero. As a result of these fluctuations, elements of 
fluid on any one of the mean lines of flow associated with the Madelung fluid 
are continually moving in an irregular way to other lines of flow, and this 
process tends to produce a more or less random « mixing » of the fluid. Like 
one of the fluid elements, the body-like inhomogeneity, which is carried along 
by the fluid, follows a very irregular path. It is then quite easily shown that 
a statistical ensemble of such systems with an arbitrary probability density, 
P(x), of inhomogeneities eventually decays into one with P=|¥Y)?. Thus, the 
irregular fluctuations in the motions of the Madelung fluid provide a model 
explaining in a natural way a possible physical origin for the statistical dis- 
tributions of the quantum theory (°). 


2. — Introduction of ‘‘Intrinsie Spin’? in Hydrodynamical Model. 


In the present paper and in a subsequent paper, we shall extend the hydro- 
dynamic model to a treatment of the Pauli equation for a spinning electron. 
Now, the Pauli equation is (1) 


whom de ea VAR E eh 
(3) th at > (v— ie.) Yin PF, ang VP, 


where Ÿ, represents a two component spinor with index, a, while A is the 
electromagnetic vector potential, V, the total external potential energy (electro- 
magnetic and otherwise), and $ the magnetic fleld. 

As in the case of the Madelung model of Schrédinger’s equation, we must 


(*) In a previous paper (see third paper of reference (*)) a similar theorem was 
proved for the model in which the W function is assumed to represent a field of force. 
In this case, the random fluctuations come from perturbations arising outside the system 
under investigation. 

(!°) We do not include the spin orbit coupling term, 6: PX€, because a consistent 
treatment of this term requires a relativistic theory, since it is of the same order of 


magnitude as the Thomas precession. This term must be dealt with in terms of a causal 
interpretation of the Dirac equation. 
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define a fluid density and a fluid velocity. Now, it is well known that the Pauli 
equation admits a conserved charge and current given (except for a factor 
of e) by 


(4-4) o = Ww, 
7 
(4-b) ji VP PIP) LA. 
2m? C 


Writing j = ov, we may then define the velocity as 


LE eee 
2mi CR DRS ARE: tre 


(4-c) vi 
It is then evidently consistent to assume that the fluid density is given by 
o=‘P*P, since the conservation equation obtained from the Pauli equation 
can then be interpreted as describing the conservation of fluid. 

o and v do not, however, permit a complete definition of the physical meaning 
of the Ÿ function; for since Y is a spinor, it contains four independent quant- 
ities, as indicated below 


a + d 
>) Sa ( + ” ; 


where. the a, b, c, d are real. On the other hand (4-a), defines only 
o= a2b?+¢?+d? while (4-b) defines only the derivatives of the a, b, e, d. 

In order to define all parts of # more directly in terms of physical pro- 
perties of the fluid, we shall therefore have to assume some new properties 
for the fluid. Now, since the Pauli equation deals with the electron spin, it 
seems natural to assume that our fluid should have a new property, which we 
may call an intrinsic angular momentum connected with the spin. By this, 
we mean that the total angular momentum density of the fluid should include, 
in addition to the «orbital » contribution of mo(r xv), an additional contri- 
bution depending on some parameters connected with the internal motions 
of each fluid element. To obtain a possible physical picture of where such an 
intrinsic angular momentum could come from, we may suppose, for the sake 
of discussion, that the fluid is constituted of molecules, and that the spin 
motion of the constituent molecules may contribute to the total angular mo- 
mentum of the system. On the other hand, the same result could be achieved 
in a much more general way. For if the fluid has an inhomogeneous structure, 
then these inhomogeneities will in general have a certain inertia. If a fluid 
element is turning, the inhomogeneities will turn with it and contribute to 
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the total angular momentum. Such inhomogeneities might, for example, be 
stable or semi-stable pulse-like structures formed in the fluid itself, or they 
might be small highly localized vortices or eddies. As far as our purposes in 
this paper are concerned, however, the origin of the intrinsic angular mo- 
mentum of the fluid is irrelevant. All that is relevant is the assumption that, 
for one reason or another, such an intrinsic angular momentum exists. In 
order to have a convenient model in terms of which we can work, however, 
we shall assume in this paper that the intrinsic angular momentum is 
due to the turning of very small quasi-rigid bodies of which the fluid 
is supposed to be constituted (just as ordinary macroscopic fluids are 
constituted of molecules). We shall see that this simplified model is adequate 
for giving a causal explanation of the Pauli equation. On the other hand, 
there is no reason inherent in the model why the bedies must always be rigid; 
or, indeed, even why the inhomogeneities must necessarily be regarded as 
arising in distinct bodies out of which the fluid is supposed to be constituted. 
Nevertheless, under the conditions in which the Pauli equation applies, we 
shall assume that the fluid acts so nearly as if it were composed of distinct 
rigid spinning bodies, that we may use this model as a simplifying abstraction 
as one, for example frequently simplifies the treatment of atoms by replacing 
them by idealized mass points even though they really have finite sizes. 


3. — Kinematie Description of Rotations. 


We shall now proceed to develop a kinematic description of the rotations 
of a body, which is particularly easy to apply to the Pauli equation, because 
it works in terms of spinors of the same kind as appear in this equation. 

The first problem is to specify the state of rotation of a body. This can 
be done in terms of the three Euler angles, 0, g and y, where 0 represents the 
angle of the principal axis (1) with a Z axis fixed in space, g represents the 
angle that the projection of this axis makes on the X-Y plane, and y represents 
the angle of rotation about the principal axis (1) relative to the intersections 
of the plane of the principal axes (2) and (3), with the X-Y plane. 

We wish now, however, to connect these angles to a spinor. To do this, 
let us imagine that we always start with the body in a standard orientation, 
in which its principal axis (1) is directed along the Z axis, which we have fixed 
in space. We then make a rotation, R(0, 9, y), which carries the body from 
its standard orientation to its actual orientation. This rotation can be carried 
out in three steps. First, we make a rotation, R,(y) through an angle y about 
the Z axis fixed in space (which is at this time also the principal axis (1) of 
the body). We then make a rotation R,(0) about the X axis, of an angle, 0. 
Then we make a rotation R;(~) again about the Z axis fixed in Space (but this 
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time the principal axis of the body is no longer parallel to the Z axis). If the 
reader will draw a diagram, he will readilly convince himself that these ro- 
tations give just those described by the Euler angles (!). 

Let us now carry out these operations mathematically. We shall assume 
that the state of the body at rest is to be represented by what we shall call 


: à ch = a 
the standard unit spinor (7?) fo = (o . When we carry out the rotations des- 
cribed above, we shall see that a rotation involving arbitrary angles, 0, 9, y, 
ean lead to an arbitrary unit spinor (such that B*B—1). Thus, any unit spinor 
ean be interpreted as a specification on the angles of rotation of a body. In 
terms of spinor notation, the first rotation R,(y) is represented by the matrix 


exp [ic,p/2] = cos y/2 + io, sin v/2 . 


Applying this to the standard unit spinor, we get 
| 4 sew Sate oye 
(6) Ry (y)Bo = exp [0/2] 0) — exp liw/21|)- 


We now apply the rotation, R, = exp [0.0/2] = cos 02 + io, sin 0/2., We 
obtain 


| fe cag AT 0/2 
(7) Rep = exp [ty/?] \ sin 0/2 COS DI RA i sin 0/2)” 


_ Applying R, = exp [toz/2)], we get 


pa _ (008 6/2 exp [i(p + )/2)\ _ i‘ 
(8) PRE (| sin 0/2 exp a) fs (a 


Writing 


Moio 
(8-a) p Le (; ale A ’ 


(1!) For a clear diagram of the Euler angles, see H. GOLDSTEIN: Classical Mechanics 
(Cambridge Mass., 1950), p. 107. | 

(12) It will become clear that the spinor that wé associate to the standard orient- 
d that all that is important is the relationship among spinors 


ation is arbitrary, an 
We have chosen the above standard spinor 


corresponding to different orientations. 
for the sake of simplicity. 
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we get 
{ b, = cos 6/2 cos (y + )/2, b, = cos 0/2 sin (y + 9)/2, 
oa | b, = — sin 6/2 sin (y + 9)/2 , b, = sin 6/2 cos (y — œ)/2 - 


The b, evidently satisfy the relation, 


(9) Bi + + B+ = 1. 


The b; are just the well-known Cayley-Klein parameters of the rotation group (1%); 
of which evidently only three are independent. By specifying the b; we can 
solve for the Euler angles, 0, g, y. Indeed, the specification is two-valued, in 
the sense that there are two values of the b; for each set of Euler angles. 

Thus, we see that a unit spinor can be interpreted as defining the orientation 
of a body in space. We shall therefore tentatively interpret the spinor appearing 
in the Pauli equation in this way, and then show that a consistent interpret- 
ation can so be obtained (1). 


4. — Canonical Relations for the Cayley-Klein Parameters. 


Now, we shall set up a classical canonical formalism, in which the Pauli 
equation is used to define the equations of motion of the Cayley-Klein para- 
meters, and therefore the equations governing the rotations of the bodies 
making up the fluid. 

Our first step is to write down a classical hamiltonian from which the 
Pauli equation can be derived. This hamiltonian is formally the same function 
of the spinor Ÿ as that which appears in the hamiltonian from which the 
Pauli equation is derived in the usual form of the quantum theory. This is 


QE im DCR i eh is 
(10) 4 ee (9 ie.) + y PES 760108 CE 


The above hamiltonian will lead to the correct equation for Ÿ provided that 
we assume the classical Poisson bracket relations 


(11) (PE (a), Vla] = 0 a(x — +). 

(**) See H. GOLDSTEIN: Chap. 4. Note that we have rotated from body to space 
axes, while GOLDSTEIN rotates from space to body axes. 

(44) Spinors have already been used to describe the orientation of a rigid body. 
See, for example, H. B. G. Casimir: Rotation of a Rigid Body in Quantum Mechanics 
(The Hague, 1931). 
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For then we obtain just the Pauli equation 


h? 


2m 


L ie 5) PI VP, + DESIO 


SOSIA TRO A 
ot 


2me 


where we have used the fact that the Poisson bracket of Y with an integral 
such as H is defined with the aid of the well-known functional derivatives. 
Thus 


SH OF 6H OF - OH OF , 6H -dF 
9 Cp a ie ao 
(12) LH, P= jp, ope SP, DPR OWE OW, OPE OW, 


In order to simplify the treatment, we now split the whole space into 
volume elements so small that Ÿ does not change appreciably within them. 
We may then define Y/(x,,) as the mean value of Ÿ in such a region, centering 
on the point x,. By integrating eq. (11) over small regions of x and x’. 
centered respectively at x, and x, and dividing by (AV), we get 


dad O(Xms Xn) 


Ee) YT (xe) | == ah, (AV) 


| We shall now find it convenient to introduce a new spinor, 


m3) y= VRAVY. 
We obtain 
(14) [xt (Xin) Av(Xn) | “= ga VOX ms Ce ae . 


We then write 


(15) ee La ie = Li : 


(a3 + 143) 


From the Poisson bracket relations, (14), we then deduce that 
(di, a] = [43, al=1, 
(16) 


[a,, as) = [@, Gs] = [@,, a5] = [42, a] = 0. 


Thus a, is the momentum canonically conjugate to 4, and a; is canonically 
conjugate to ds. The a, are evidently proportional to the Cayley-Klein para- 
meters, b, (given by ed. (8-a)). Thus, we have found the canonical relations 
among the a;, implied by the Pauli equation. 
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5. — Definition of Spin Angular Momentum. 


In order to motivate further steps in the interpretation of the Pauli equat- 
ion, let us recall that in the usual quantum mechanics, the total spin is given 
by (#/2) [ro dx. This suggests that (h/2)(¥*oW) is a spin density, and 
that S = (h/2)(V*eP)AV is the total spin in the small element of volume AV. 
Expressing Y in terms of 7 through eq. (13), we get 


(17) S = }(x*ox). 
The total number of particles in this region is just 


(18) AN = W*WAV= 7*z/h. 


Thus, the spin per particle is 


RV SOY IO 
2 y%y 2 PP 


(19) S— 


We see then that the maximum spin per body in a given direction is #/2. 
This is in agreement with what one gets for the spin « observable » in the usual 
quantum theory (15). 

With the aid of the P.B. relations (11), it can be shown be means of a 
simple calculation that the components of the spin vector, S, satisfy the cyclical 
relations, characteristic of angular momenta 


(20) [Se S,] = 8,. 


Thus, the Poisson-Bracket relations obtained in the derivation of the Pauli 
equations from a hamiltonian are just what are needed to lead to the correct 
P.B. relations for the angular momenta defined in eq. (17). 


(13) In our theory, the spin vector is a continuous variable, with arbitrary projections 
on any axis, and with a total magnitude of s2—#%2/4 (as can easily be proved by using 
the Pauli identities). On the other hand, the spin « observables » appearing in the 
usual quantum theory may have only certain discrete projections on any divo axis 
while the magnitude of the «observable» for the total spin is #42. The conn of +. 
difference arises in the circumstance that in the model that we are propos here 
the spin « observable » is, as we shall see later, a kind of statistical or over-all propert 
of the motions of the bodies in the fluid. Thus, there is no reason why it should È 
identical with the angular momentum of one of the bodies in the fluid, although there 
will in general of course exist a certain connection between the « observable » and the 
spins of the bodies, which we shall discuss in more detail in the subsequent paper 
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In the present theory, the basic spinor, x = 2 È Si contains only 
V2 \a3 + ia, 

two canonically independent variables. We shall now see the physical meaning 

of this reduction of the number of basic variables to two. To do this, we find 

it convenient to go to a new set of canonically independent variables which 


we take to be two of the Euler angles, y and gy. In terms of the a,, these 
are (as can easily be seen from equ. (8-b)) 


Ho 0 - 
(21) y — te! nae + tg} = : p E to! SER le Da 


By means of a simple calculation, we prove that the quantity canonically 
conjugate to — y/2 is oh, which is proportional to the total number of part- 
icles in the region (see equ. (18)), while the quantity canonically conjugate 
to —w/2 is S,, which is equal to the Z component of the total angular mo- 
mentum of these particles. Thus, we have 


(22-a) no-5 = (et ay + as “ts a) (te (UP) + ts “ — À | ; 


2 2 ay dy 
2 2 2 2 
La (ai À da dg da) ol da ol da == 
Sai (se aprano; » (te di ve dal] = 


But now we can show by a simple calculation that is also equal to the 
angle g' made by the projection of the angular momentum vector S on the 
X-Y plane. Thus 


; È ID ol: de BETA 
ì ; x 
Pro, ds — agli 1 0/\03 + 004) _ Wola dits, 


ee a Pro E a cat ale Te al ap 4,0 + Agha 


da — t0,)]\1 0/\% + ia, 


Now, from eq. (21) 


go = (2 2) (a 428) = 
Di Wy 4} \ Ada MA + Ada 


Thus, we see that g'= y. Moreover, the co-latitude angle, 6’, of the angular 


2 2 2 2 

S a, + Ge — dg à 

: I 4 1 2 3 4 B { fr 

mentum vector, is defined by cos piede A ut from 
Ru i S ay + a + as + da 


equ. (8-b) we can calculate the angle 0, made by the principal axis (1) of the 
i 7 axis. This is given by cos? = de aie rds = a Thus 
body with the 7, axis. This 18 giv PIE tate git gt =e ; 


DZ 6 
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We conclude then that in the Pauli theory, the angular momentum, S, 
is always pointing along the principal axis (1) of the body. Such a special 
orientation of the angular momentum cannot be maintained for the most 
general kind of torque that may act on the body. In the next section and in 
a subsequent paper, we shall see, however, that the Pauli equation implies a 
special kind of « quantum-mechanical torque » that permits this condition to 
be maintained as a consistent subsidiary condition. 

The special relation between the direction of S and the principal axis (1) 
is what permits us to reduce the number of independent canonical variables 
in the theory. When this relation exists, our theory becomes equivalent to 
a theory of the angular momentum of a point dipole, which has already been 
treated by several other authors (**"*). In our subsequent paper, we shall 
see, however, that the theory developed here can be generalized to treat cases 
in which $ does not necessarily point along the principal axis. 


6. — Relationship Between Spin and Velocity. 


Let us now recall that equ. (4-c) defines the velocity of a particle in terms 
of our spinor. Since the spinor is already interpreted in terms of the orientation 
of a spinning body, eq. (4-c) clearly implies a certain relationship between 
the spin and the velocity. We shall now study this relationship and show 
that it can be understood in terms of a reasonable physical model. 

We begin by expressing our spinor (with the aid of eq. (7)) as 
= a he SS ly da poi 5 where È is real. 
à sin 0/2 exp [i(y —o)/2] 


Using this value of Ÿ in eq. (4-c), we readilly get 


h e 
(23) DES A A li IN 


The above, is however, an expression of the velocity in a form that has 
already long been familiar in classical hydrodynamics. Indeed, in classical 
hydrodynamies (!°), it is shown that an arbitrary velocity field can always 


(8) G. F. UHLENBECx and S. GOoUDSMITH: Nature, 117, 264 (1926). 
(2) H. A. KRAMERS: Quantentheorie des Elektrons, (Leipzig, 1938). 
| (5) C. M. LATTES, M. ScHENBERG and W. SCHUTZER: Anais da Academia Brasi- 
leira de Ciéncias, XIX, no. 3 (September 1947). 


(1°) See, for example, H. Lams: Hydrodynamics (Cambridge, 1953), p. 248. 
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be expressed as 
(24) Dek yy A 
C 


where È and 7 are scalars, called the Clebsch parameters, and S is also a scalar. 
The curl of the velocity is then 


(25) Vxv=VexW—-VxA. 


(This is quite a different expression from the more common form 
v = VS'+ Vx B— (e/c)A, in the sense that S’ and S are different functions, 
‘since div (£ Vn) is not in general equal to zero). 
The Clebsch parameters have recently been used by Dirac (?°) and 
and others (21) to treat a theory of a classical electrified fluid ether. They 
‘have also been used by TAKABAYASI (*) to treat possible vortices in the Ma- 
 delung fluid, and by SCHONBERG (7) to treat vortex motions in a more general 
context. 
In order to establish a basis of comparison with the Pauli theory, we shall 
first treat the ordinary classical hydrodynamics of a charged fluid, capable 
of maintaining a pressure gradient, in terms of the Clebsch parameters. One 
of the advantages of these parameters is that they permit the formulation of 
hydrodynamics in terms of a variational principle. Indeed, if @ is the density 
of the fluid, then the hamiltonian is just equal to the total kinetic plus po- 
tential energy of the fluid (including the potential energy due to compression) 


(26) H.= ig 


£ (vs LéEVy— <A) + oem + ID da, 
where A is the vector potential, y, the scalar potential, e/m the ratio of charge 
to mass for an element of the fluid in question, and f(g) defines the pressure 
| through P= 0f/00. 

Now if we adopt essentally the same Poisson Brackets as were derived 
from the Pauli theory; i.e. 


(27) [o, S] = dla — x"); [(o€), n) == x — x"); 


(2) P. A. M. Dirac: Proc. Roy. Soc., A 209, 291 (1951); 212, 330 (1952); 223, 


439 (1954). 
(21) See, for example, O. BUNEMAN: Proc. Camb. Phil. Soc., 50, 77 (1954). 
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se 


then we obtain the following equations of motion 


(28-a) > + div @ (vs RAI 24) = 2 + div (0v)=0, 
(28-b) 28 glow + St At tot i Lo, 
(28-c) = = 

(28-d) si + (0-Vin = di = 0 


Combining (28-b) and (28-d) we get, 


08 | e (VS + €Vn— (e/e)A)? | <> fil o 


GU 
Ge) BES CT ee 2m m Co 


Eq. (28-a) is just the equation of conservation of fluid. Eq. (28-b) (and 
therefore also (28-e)) is a generalization of the equation for the velocity po- 
tential, S, which enables us to take vorticity into account (through & and 7), 
Eqs. (28-c) and (28-d) are very interesting, for they say that if we follow a 
moving fluid element, £ and 7 are constants of the motion (but not at a point 
fixed in space). Thus, the surfaces, é(a, y, 2, t) = const. and 7 = 7(@, y, 2, t) = 
= const., define the tubes of flow of the fluid. 

To understand the motion in more detail, let us focus our attention on 
eq. (28-e). This is analogous to a Hamilton-Jacobi equation, but the equation 
involves, in addition to the electromagnetic potentials, and the pressure po- 
tential (P = 0f/00), a set of equivalent potentials that are functions of the 
Clebsch parameters (??). In fact, the Clebsch parameters contribute an ad- 
dition of ?(07/0t) to the scalar potential and — éVy7 to the vector po- 
tential. We may thus conclude that formally they contribute to the equation 
of motion of a fluid element a « pseudo-Lorentz force », given by 


(29) RO: 


where the « pseudo-electric » and « pseudo-magnetic » fields are given 


\ a 


Poe on 0 0Ë on 
c= v(e3) L ape Nava; 


§'= —Vx(EVn) =— VEXVn. 


(22) We are using here an argument due originally to M. ScHòNBERG: Nuovo Ci- 
mento, 11, 674 (1954). 
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Thus, for the additional force, we get 


DE 
(30) F'= 2 Vy 
ct 


4 VE vx(VExn) - 


| © 


Writing — v x (VEX Vn) = (v: VE) Vn— (w:Vn) VE, we obtain 


MICA en dé dy 
PF = — + È — — . È È a 
3 + (v-V) a (@-Vig= qi Vn at VE 


Since from (28-c) and (28-d) we have dé/dt = 0 and dy/dt = 0, we get HQ, 
Thus the Clebsch parameters formally add to the vector potential a set of 
quantities that produce a pseudo-electromagnetic field for which the pseudo- 
Lorentz-force vanishes. Therefore the Clebsch parameters appearing in 
eq. (28-e) do not alter the equations of motion of a fluid element (23). What 
they do is to enable us to consider in a canonical formalism an ensemble of 
trajectories for the fluid elements which have a vorticity that does not come 
from the effects of the electromagnetic vector potential, A. Thus, they make 
possible a generalization of the Hamilton-Jacobi theory; for in the latter, 
we consider only ensembles in which mv — (e/c)A is derivable from a potential 
(equal to the action function). 

Let us now return to the Pauli theory. Our first problem will be to 
compare the classical hamiltonian (26) with the Pauli hamiltonian (10). By 
means of a simple calculation, we obtain the result that the Pauli hamiltonian 
can be written as the sum of the following three terms: 


H= Hy + Hr Ho 
where we have 


he ,[Vy Vaenren\ 
(31-4) farsi rs | 5) Co80 5 | Axe 
h? (Vo')? 
- H, = |— > 4 
(31-5) È (E 0 dep 
( ne i È 
(31-c) He | dx E o'((V0)* + sin? 0(Vp)?) + 
+ a (cos 09, + cos 0 cos pa + COS Asin pDy\ , 
2me 


where § is the magnetic field. 
(23) Note that the fact that dS/dt=dn/dt=0 plays an essential role in making 


this result possible. Later, in connection with spin theory, we are going to have 
e derivatives; and then, the Clebsch parameters will imply 


non-zero values for thes 
additional terms in the force on a particle. 


5 — Supplemento al Nuovo Cimento. 
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In the above equations, we have, as in eq. (4), defined = PRO; 
where o' is the density of bodies. 

In order to compare the Pauli hamiltonian to the simple hamiltonian, He, 
of classical hydrodynamics, we note that 9’ is canonically conjugate to —y/2, 
(o’ cos 0) to —q/2. Thus, if we replace @ by o’, and write 


(32) E=così, nih=— 9/2,  Sfh=— [2 


then H, becomes equal to kinetic energy term in eq. (26). 

The terms that are left over in the Pauli equation are then H, and H,,. 
We shall see in eq. (34-b) that H, leads to the « quantum-potential » term of 
the Madelung fluid. Thus, the two terms, H,+H, represent the energy of 
a charged Madelung fluid, in which vortex flow as well as potential flow are 
taking place. 

As for the third term, H,, , it is clear first of all that it contains a part 
describing the energy of a spinning dipole in a magnetic fleld $. The magnetic 
moment of this dipole is eh/2me. 

If we assume that our bodies have a radius of the order of the classical 
electronic radius, 77>=2.8-10-" em, then the magnetic moment is u=Aergofe 
where v is the velocity of rotation of the body at its periphery, and where À 
is a constant of the order of unity depending on the charge distribution in 
the body. We then obtain 


v h 137 

c Amor CAS 
Thus, the body must, at its periphery, be moving much faster than light 
In a non-relativistic theory, such as we are discussing in the present paper, 
this, of course leads to no difficulties. In a later paper, concerned with the 
extension of the theory to the Dirac equation, we shall see that the spin 
motion becomes closely coupled with the mass motion of the fluid, so that a 
spin anguiar momentum of eh/2me for the body is possible, even when the 
speed of its periphery does not exceed that of light, because the fluid motions 
are further connected with the body rotations by a spin-orbit coupling, not 
present in the non-relativistic theory. 

As for the remaining part of H,,, this represents an interaction between 
the directions of the spins of neighboring bodies, which functions essentially 
as a spin-dependent addition to the « quantum potential ». To obtain a more 
definite model for this part of the spin energy, we may let dA represent the 
total angle between two spin vectors, separated by a distance, dx. Then, by 
going to spherical polar coordinates in the space of the spin vectors, we can 
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show that 


\dA 


3=1V0-dx + sin? 9|Vp-dx|?. 


Thus, the energy, H,, is just an average of (6A), taken over a small region 
that surrounds the body and weighted with the density, o. We can therefore 
interpret H,, physically as the result of an assumed short-range interaction 
between bodies, which tends to line up the directions of their spins. For ~ 
example, the energy of two neighboring bodies could be proportional to 


(33) Uinr, = — C08 (0A) = —1 + (04)°/2. 


Let us now obtain the equations of motion that follow from our hamil- 
tonian. To do this, it will be convenient to write H as 


oe e) he (Vo) 
a= | Li (vs LEVy “Al dx 4 mI o dx 4 | emsax, 


| where, 
| h? : € he AVE), 3 0 La 
GE, = > [(V0)*+ sin? 0(V@)?] + CA syY= Di Soa) (1 SU ss me © . 


The equations of motion then become 
| 


| Co! , 
Meta) 3 divov—0, 


ot 
| os me uo (Vo) 
| (34-b) ii sa E(v-V)n | 9 malo 9 (o') | Hs 0, 
34 dee 
| Ee?) fea De OT) 
(34-d) rta 


Eq. (34-a) is then just the conservation equation while eq. (34-b) is the 
pseudo Hamilton-Jacobi equation for a fluid with the usual quantum potential, 
! 1\2 i 
a di ve ze wer) and a spin-dependent quantum potential, H,— 
2m 0’ 4m (0°)? 
a Ay Vin: . . . . 
Eqs. (34-c) and (34-d), which tell us how the spin directions change with 
time are in precisely the same form as that of the equations expressing the 


64 D. BOHM, R. SCHILLER and J. TIOMNO 


torque acting on a point dipole, derived previously by SCHONBERG (*) in a 
purely classical problem. 

We may now, however, obtain another instructive form of the equations 
for the spin variables by noting that H,, can be written in the form 


he È Si ces 
(99,0) nio 2m 2 ca amo” 2. 
È 2 | S| Ta i 4 
Noting that 0 = a we also obtain the expression 
“Thi 7/08 i ME 4 
3 5 = | |— — 1 — (S: da. 
(39-6) H+ Hu = [| 5) /e +5 (5-8) 


We can obtain the equations of motion for s, by writing x = S/o, and 
using the P.B. relations for the angular momentum. We get (with the aid 
of eq. (35-a)) 


ds s lo os e as 
Hp —— — — x). 
=o) dt 0 S > 0%; (e =| a ME (839) 


Thus, the spin vector precesses about the magnetic field, with angular fre- 
quency, eS/me, and there is an additional precession with angular velocity 


eee Os | 
Dr (o =), which results from the torques produced by the neigh- — 
i ON; i 


boring spins. In some ways, one may think of this extra torque as due to -" 
a kind of « quantum-mechanical » addition to the magnetic field. But to make 
the analogy correct, one must think of the magnetic field in a polarizable 
medium where each dipole interacts strongly with its neighbors. 

Finally, we shall discuss the motion of an element of fluid. This is best 
done in terms of the expression for the energy-momentum-stress tensor. To 
derive this tensor, we start with the lagrangian leaving out the vector and 
Scalar potentials, which do not alter the problem in any essential way 


£ "ras Monza he (Vo)? 
3 = RE | DR == 
(37) L | 0 È mica 8 (Vy EV) Sn (0)? 
h? 
— a ((V6)?+ sin? 6(V@)*)| dx. 


It can easily be shown that this Lagrangian leads to the correct canonical 


relations (27) between € and 7, and between 0’ and S, and that it leads to the 
Pauli hamiltonian (10). 
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Now the canonical energy-momentum stress tensor is given by (?4) 


hs ONCE 


OL 
Ti; = AIA$ IA ; Toi = = 
2 O( Of,/Ox;) Ox A 2 O( Of, Ot) Oa, ? 


3 
LE > A(Of,/ot) of ’ 


where the f, represent all the possible field quantities. Ti; represents the 
momentum density, and 7» the energy density. We readily obtain 


oS 


lo) 
To =e (+65), 
da 21 


reel e (È je cd) LOSCO 


A Î 
Am 0%; 0%; 


# 2 Dy 
i JE 0 Limo ÎE BLY 


/ n a 7 ey 
4m OX; 0%; Di 0%; 


| where È and 7 are given by eq. (32). 
Pe IL oS 01 
The energy per particle is then 5 + Le , which agrees with what we 


obtain with the « Hamilton-Jacobi » equation (28-e). The momentum per 
| particle is just p= VS + éVy, in agreement with what we have been as- 
| suming. The stress tensor, T;;, contains three parts. The first part is just 
mo'v;%;, which is just the usual term representing the effects of mass motion. 
The second part corresponds to the quantum potential and the third part to 
the effects of the interactions between the spins of neighboring particles. The 
equation of conservation of momentum can then be written as 


D 


ot; 


- (0/0: 0e 
(g's) +2 a, 


ind) 


With the aid of the equation for conservation of particles, and the relation 


We É 00" ce) 0 (6) no A 


do 0%; © On; dx; A ox of dx; ’ 
we obtain 
dp ov dA il a x 
(39) ‘ 7a m (= + (D Vie) —— VU : (08;;) ’ 


(24) See, for example, G. WENTZEL: Quantum Theory of Fields (New York, 1949). 
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We see then that the fluid element moves under the action of the « quantum- 

h> do do 
4mo' 0x; dx; 
and an additional quantum- ont contribution to the stress tensor, arising 
from the interaction of neighboring spins. 


mechanical » stress tensor, , which leads to the quantum potential, 
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1. — Introduction. 


In a previous paper (!), to which we shall hereafter refer as (A), we have 
_ developed a causal interpretation of the Pauli equation, in terms of the model 
| of a fluid composed of spinning bodies, in which the spin angular momentum 
of each body is directed along its first principal axis (i.e., its axis of symmetry). 
The object of the present paper is two-fold, first, to generalize the above model 
to show how the Pauli equation results as a special case of a theory in which 
the angular momentum may point in an arbitrary direction, and secondly, 
to illustrate our model of the Pauli theory in terms of the simple example 
of a stationary state of an electron in an atom, and to apply this example 
in the theory of measurements. | 
In Secs. 2 and 3 we shall develop our more general model and show that 
in terms of it, the Pauli equation follows as a consistent subsidiary condition. 
In other words, if the angular momentum points along the first principal axis 
at any instant, say t= 0, then the torques will be such that this condition 
will be maintained for all time. Moreover, small deviations from this con- 
dition lead to a rapid oscillation of the motion about a mean, in which the 


(1) D. Boum, R. ScHILLER and R. Tromno: Nuovo Cimento, 1, 48 (1955). 
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condition is satisfied. Physically, such an oscillatory motion corresponds to 
a precession of the angular momentum about the direction of the axis of 
symmetry, in which the component of the angular momentum perpendicular 
to this axis averages out to zero. Thus, even in the case of a general orientation 
of the angular momentum vector, the Pauli equation would still be a good 
approximation for processes that are slow compared with the rate of pre- 
cession described above, while for more rapid processes, it would cease to be 
a good approximation. Instead, a more general equation (which will be seen 
to be non-linear) would have to be used. Thus, we are led to a specific example 
in which the causal interpretation of the quantum theory implies new kinds 
of equations in connection with high frequencies (and therefore high energies), 
equations which would reduce approximately to those of the usual quantum 
theory only at those relatively low energies for which the usual theory is 
known to be valid. 

In Sec. 4, we then illustrate our model of the Pauli equation in terms of 
a stationary state of an electron in an atom. This illustration serves not only 
to bring out important topological properties of a field of body orientations, 
such as is implied by our model, but it also shows how a continuous distri- 
bution of spin angular momentum throughout the fluid can lead to discrete 
or « quantized » possible values for the angular momentum of the fluid as a 
whole. Quantization of angular momentum (and of other variables) is there- 
fore seen to arise as an over-all property of the fluid, on the basis of a lower 
level of continuous motion, much as discrete frequencies of stationary modes 
of vibration arise in connection with the motions of continuous fields in 
classical mechanics. 

In Sec. 5, we indicate how the process of measurement is to be treated in 
terms of our model. What is observed in a measurement of the angular mo- 
mentum «observable » will be seen to be just one of the discrete possible 
stationary values for the total angular momentum of the system, and not the 
continuously distributed angular momenta existing in the various parts of 
the fluid. 

Finally, in Sec. 6, we summarize the essential features of the model, and 
suggest possible further directions of research. 


2. — Treatment of Arbitrary Direction of Spin Motion for a Fixed Rigid Body. 


Our first step will be to extend to the case of arbitrary spin motions the 
hamiltonian formalism in terms of spinor variables developed in paper (A). 

To do this, we begin with the case of a rigid body at rest. In the next 
section, we shall then extend the theory to the more general case of a field 
of rigid spinning bodies with translational motions. 
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The lagrangian of a rigid body with no translational motion is just the 
kinetic energy due to its spin, which is (?) 


(1) a 1 5 TOW, , 


ik 


where the ©; represent the components of the angular velocity, and /;, is 
the tensor for the moment of inertia of the body. The components of the 
angular momentum are then given by 


(2) Sage DI VESTOT 
x 

so that 

(3) n= LD sw. 


The next problem is to express the ©; in terms of spinors. Now, the 
relationship between the ©; and the Euler angles is (*) 


D, = ÿCOBD= EP, 
(4) ©, = Ÿ Sin 0 sin y + § cos y, 


i), = y sin 0 cos y — Osiny. 


We now return to our spinor f = (al defined in paper (A), eq. (7). We have 
2 


3 (wp + AR NERO SATO 
Bu exp fe) (5 (y + ®) cos Bags sin 3) 


: IZ ATO 8) 0 
Ba = è exp BS (Se P) sin 5 + 50085) 


Now let us consider the quantity 
7 È, È, ald . . . NO, + PISA >) = i(p così + D) 
Bro. Bro,B=i (cos? à (p +o) — sin? 5 (p— d))= il ). 


Hence we obtain 


(i= i(B*0.f — p*o:P) 


(2) In the present paper, we leave out the electromagnetic field, which makes no 


i i i heory. 
essential change in the formulation of the t Le: o 
() H. GOLDSTEIN: Classical Mechanics (Cambridge Mass., 1950). See Chap. 4. 
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And more generally, because of rotational symmetry 
(5) w = i(b*of — f*of). 


The kinetic energy then takes the form 
A doo 5 N é 
(6) LE a) — p 4o4s) 6) 


We now define the momenta canonically conjugate to the components of 
the spinors f and £*, which are 
OL OL : CILE OL 
ep,’ ops ape ope 


= 
ten | 
— 


Pai: 


We also define the spinors 


ila a Du 
(8-4) T= (i) A ta , 


so that formally speaking 
(8-b) Me = * — 


We then obtain from (6) 
(9-a) a =—t(0's)p*s 
(9-b) [a —— KC RRS one 


Note that in eq. (6), s is a function of È given by eqs. (2) and (5). As can 
easilly be verified, however, this dependence merely leads to a multiplication 
by a factor of two of the result that would be obtained if s were not a 
function of f. 

We readilly see that under rotations, 7 transforms contragrediently to b. 
To prove this, we form the quantity x. From eq. (9-a), we obtain 


af = — i(B*(0-s)f), 


which is evidently a scalar, so that x must transform in the same way as pr, 
or in other words, contragrediently to ft. 
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To solve for s;, we multiply (9-a) by fo;, and (9-b) by p*o;. We obtain 


i 
(10) se) (108 — x*oB*). 


We shall now find it convenient to obtain the expressions for the projections 
of the angular momentum vector on the principal axes of the body. It is evident 
that each of these projections is a scalar with regard to rotations of the space 
axes, since by definition, these components will not depend on which space 
axes we choose for the expression of s. We can therefore evaluate the com- 
ponents of the angular momentum along the body axes by choosing our space- 
axes to agree with the body axes. In this case, according to paper (A), eq. (7), 
we have p= (3) and for the components of the angular momenta, we get 


— 


(11-a) TL M a EA RTE 


We also introduce for convenience, the quantity, T5, which as we shall 
see presently is zero in our problem, 


NV 7% 
(11-b) Tees È det 


We wish now to express the 7, in terms of the spinor quantities, taken 
in an arbitrary frame of reference. To do this, we seek a set of scalar functions 
of the spinors x and f which reduce to the 7, when we chose the space frame 
to be the same as the body frame. 1 

Let us first consider the scalar 78. In the frame in which fp = Ge this 


reduces to 


np =m, = T_T 


so that 
Qk Ook 
(12-a) T2 (rp ka ) 
nb + a*p* 
(12-b) T, = ues “ie | 


Now from eq. (9-b), we see (since s is real) that 7, = 0. This relationship 
is a subsidiary condition whose meaning can be seen by writing f and z in 
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extended form 
PR ti air De) bag — ei sa 
PE ep i 


where p; is canonically conjugate to b;. We then get 


(1)3 mp + a*p* = bp, + Dep. + bss 4 baPa + 


If we consider the space of the Cayley-Klein parameters introduced in 
paper (A), eq. (8-b), we see that from the relation bî+05+b3+0f=1, it follows 
that we must remain on a unit hyper-sphere in this four dimensional space. 
Now eq. (13) then expresses the fact that the «radial momentum » in this 
hyper-space is zero. Such a result is to be expected, since the kinetic energy 
is a function only of the angular velocities, which do not involve the « radial » 
component of the velocity in the C.K. space, so that the corresponding mo- 
mentum must be zero. 

From (11-a) and (11-b), we readilly obtain for the total spin 


(14) S2 = T*= (ata) T° = (n*x)? (since T, = 0). 


0 


We now wish to obtain expressions corresponding to (12-a) and (12-b) for 
T, and T;. To do this, we first consider the spinor 


Û (A) 


The spinor È has the same transformation properties under rotation as 
does f itself, as can be readilly verified by direct computation with the trans- 
formation matrix exp [i(R-o)/2]. Indeed, we have p= i0,P*, which is just 
the non-relativistic form of the so-called charge conjugate spinor (*) which, 
as is well-known, has the same transformation properties as the spinor b itself. 

We now consider the invariant nf . If we evaluate this in a frame in which 
the space and body axes are the same, so that B = (1). then we obtain 


(using eq. (11-a)) 


(16-a) np = % = T,—iT,, 


(4) W. PAULI: Rev. Mod. Phys., 13, 203 (1941). 
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or 


9 a 7 Qk di % D % 
= zh De p : Te a= ; CB a DI a 


È a 


(16-b) Sb 


It is readilly verified by direct computation that (T,, T., Ts) satisfy the 
characteristic P.B. relationships of angular momenta, (DID) = Te, cer 
One also readily verifies that the P.B. of any component of 7, with any com- 
ponent of S are zero, as they have to be, since S is the generator of the infini- 
tesimal rotation of the space axes and since the T, are scalars with respect 
to such rotations. 

We shall now go to the hamiltonian formalism. The hamiltonian is 


(17) l= > SW, Va GL ===> >: ie — L . 


We must now eliminate the ©, from the hamiltonian. This is most con- 
veniently done by expressing the angular momentum in terms of its projection 
T, on the principal axes of the body. This gives 


TI, 
DITE QUE: ? 


(18) HE 


where I,, I., 13, are the moments of inertia relative to the principal axes. 
For the case of a symmetrical top, which interests us here, I, = I3= I. Let 
us write 1,—1/e. Eq. (18) then becomes 


I 


1 il 
Ci ee ee a he : 


2] 2 Je (hear tenet) Ta) 


(19) H 


To obtain the equations of motion for S, we note that the P.B.’s of S with 
the 7; are zero. Thus the components of the angular momentum, taken with 
respect to axes fixed in space, are constants of the motion. As for the T,, 
we note that 7, is also a constant of the motion. For Hand, We readilly 
obtain the equation of motion 


d ACE 1) 


dt (Tr Ves) ==} I T,(T3— T2) Py 


(20) 


Thus, the component of the angular momentum normal to the first prine- 
ipal axis rotates with angular velocity, © = ((e— 1)/I)T,. This means that 
the first principal axis rotates around the direction of the vector, s, with the 
angular velocity ©. Thus, on the average the first principal axis points in 
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the direction of s; and insofar as phenomena that do not change much in the 
time t= 1/m are concerned, the body acts as if its angular momentum were 
parallel to this principal axis. 

We note that when e< 1 (disc-shaped object), the energy is a minimum 
(for a given total angular momentum |Z) when 7) = T3; Le whens the 
principal axis points along the direction of s. Thus, for this case, the motion 
in which the principal axis is parallel to s will be stable. If «> 1, (cigar- 
shaped object), such motion becomes unstable. 

Finally, we shall express in terms of spinors the condition that the motion 
reduces to the type implied by the Pauli equation, in which the spin points 
along the principal axis. If this condition is satisfied, we have T, — T; = 0; 
or ap = 0. This leads to 


7, Bx = LAS ; or Ty | Ts, + Br /Br ’ 


so that 
(21) MA al la 


where k is a constant, and where we have chosen the factor of 7 to simplify 
subsequent expressions. But in order that 7, = (aB-+2*f*)/2 = 0, k? must 
be a real number. If this is the case, however, then we see that i8* is canon- 
ically conjugate to 8, which is what is essentially the same relation as 
paper (A) eq. (14), which holds in the Pauli theory. To demonstrate in more 
detail the relationship between x, f, and the Pauli spinor, we make the fol- 
lowing canonical transformation to the spinors A and B 


[ A= (kB + inx*/k)1/v2 , 
[PBS (k= tse) NE 
From the above, we readilly prove the Poisson-bracket relations 
(23) [A; B] == [A*, 24120); [PAY ANR Enr 
To satisfy eq. (21), we choose B = 0. Then the spinor A* has the same 
P.B.’s with A as does the Pauli spinor (see paper (A), eq. (14)). Moreover, 
we also have A= kf, so that f is proportional to the Pauli spinor. To express 


the angular momentum, we eliminate x and B from eq. (10), obtaining 
a 


(24) s = k(A*oA). 
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Thus, the theory is able to treat an arbitrary value of the spin angular mo- 
mentum. To obtain the same value as in the Pauli theory we must set k=h/2. 
In this way, we complete our demonstration showing how the canonical for- 
malism for an arbitrary direction of spin reduces to that of the Pauli theory 
when the spin points along the principal axis. 


3. — Extension to Rigid Bodies with Translational Motion. 


In the previous section, we developed a hamiltonian formalism that permits 
a treatment in terms of spinors of arbitrary motions of a rigid body fixed in 
space. We shall now extend this formalism to permit the treatment of a field 
of such rigid bodies, undergoing translation through space. 

To do this, we shall first write the kinetic energy of translation of the 
bodies in terms of Clebsch parameters as we did in the case of the Pauli theory. 
Now, in the Pauli theory, we found that one rather naturally obtained the 
expression given in paper (A), eq. (25) for the velocity, v = #i/2(Vy+cos 0 Vp), 
which involved only the pair of Clebsch parameters, cos 0 and g/2, along with 
the velocity potential, y/2. It is possible, however, to introduce aS many 
additional pairs of Clebsch parameters as we please in the definition of the 
velocity. Thus we may write 


(25) v = Vi + > EVN: è 


Now, it is true that the most general velocity field can in principle always 
be expressed without these additional pairs of Clebsch parameters. Neverthe- 
less, the physical conditions of the problem may often make it convenient to 
introduce such additional pairs. Indeed, as we shall see presently, the treatment 
of the problem of the motion of a body with a general orientation of its angular 
momentum relative to its principal axes is an example of just such a problem. 

In terms of the expression (25) for the velocity, the kinetic energy of trans- 
lation of the bodies then becomes 


(26) TES (VA + XE Vi) da. 


Now, introducing the P.B. relations 


(27) Lo(x), A] = [(o(x)élx)), n(x’) = dla — #7) 
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with all other P.B.’s zero, we get for the equations of motion 


00 7 
(28-a) cn + div (ov) =0, 

ys Yi GE Lao 
(28-b) ua > E,(v-V) a "o D 0, 
DORs = —— — (0). 
ERO), dt dt 


Eq. (28-a) expresses just the conservation of bodies, while (28-5) is an 
extension of the Hamilton-Jacobi equation. Eq. (28-c) expresses the constancy 
of the Clebsch parameters as we follow a moving body. It is this property 
of the Clebsch parameters that is most interesting to us here; for eq. (28-c) 
means that the &; and the 7; may represent internal properties of bodies (such 
as angle or spin variables). While these variables remain constant in eq. (28-c), 
if we add terms to the hamiltonian involving the €; and the 7;, then the 
equation that takes the place of (28-c) will tell how these properties change 
as we follow the moving body. 

The Clebsch parameters therefore provide a natural canonical treatment 
of the motions of a field of parameters that represent internal properties of 
moving bodies. 

Now, in,the Pauli theory, the current vector was expressed as 


J = 0 VÀ + (06) Ve. 


with 
£ = cos 0/2, n=— 9/2, À = — y/2. 


Let us recall that he was canonically conjugate to — y/2, heé to —@q/2. 
Thus writing he/2=J,, he cos6/2=J,, we have (remembering that J, is 
density of the component of the angular momentum in the direction of the 


principal axis, and J, the density of the Z component of the angular mo- 
mentum) 


J => TE; Vy = a Vy Si 


But now we are going to give up the restriction that Sis directed along the 
principal axis. To complete the description of the angular momentum, we 
shall introduce, J,, the momentum canonically conjugate to the Euler angle, 0. 
This is equal to the density of angular momentum about an axis perpendicular 
to the plane containing the Z axis and the principal axis (1) of the body. 
In order that the P.B.’s of the kinetic energy of translation with J, and 0 
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NI 


~J 


shall give dJ,/dt and d0/dt respectively, it is clearly necessary, however, that 
j shall be the same functions of J, and 0 as it is of J, and y, and of J, and y. 
This, we write 


j=d,Vy4+d,Vo + J, V0. 


Now, in order to obtain a simple model that approaches the Pauli case, 
we set the density of bodies equal to 2J,/h, which is proportional to the 
density of the component of the angular momentum along the first principal 
axis. This implies that all bodies have the same value of this component of 
the angular momentum. We shall see presently that this assumption is con- 
sistent with the equations of motion that we are going to adopt. We then get 


BY jh ab: jp 
(27) eee Vs de ombre 


The next step is to add the kinetic energy of spin to the hamiltonian. 
To do this, we introduce the quantities Q; = e7;, which are the densities 
of the components of the angular momenta taken along the principal axes. 
From the results of the previous section, this is 


(28) Hy = | sp (O° (e— DO ax 
(note that by definition, Q, = J,,). 

Finally, we must add an appropriate generalization of the quantum 
mechanical part of the energy («quantum potential» plus spin potential), 
which, for the Pauli case is expressed in paper (A), eq. (35-0). We write for 
the term 


2 ig SN 2 
(29) Ho d Di da 


5% 2m } ‘i 
The complete hamiltonian is then 


(30) H=H,+H+H,. 


Let us now obtain some of the equations of motion. First of all, since @, 
has zero P.B.’s with S;, we get for Qi =J,, 


2Q 


MONT) 0% 
ot 


(31) 
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This is just a conservation equation, for the component of the angular 
momentum along the first principal axis. Since the density of bodies, o, 
satisfies the same equation, we see that if we choose hoe/2=Q, at any 
instant, say t= 0, they will remain equal for all time. But ho/2 = Q, 
implies that all of the bodies have the same value of the component of the 
angular momentum along the principal axis. Since this component is a constant 
of the motion, the condition that it has the same value for all the bodies will 
evidently likewise be maintained for all time, if it is satisfied at any time. 

As for Q, and @:, we get 


d [Qs— iQ) _ (€—1)Q1 d— 102 
na 

Thus, as in the case of the freely rotating body, the component of the 
angular momentum perpendicular to the first principal axis rotates with ang- 
ular velocity @ = ((e—1)Q;)/Zo. We may evidently then adopt Q.=Q; = 0 
as a consistent subsidiary condition, so that if the body is set with S ori- 
ginally in the direction of its first principal axis, then the torques will be such 
that this condition is maintained for all time. Moreover, if the principal axis 
is nearly in the direction of the angular momentum, then @, and Qs; will rotate 
with angular velocity, w,so that for phenomena involving characteristic times 7, 
which are much less than 1/w, the body will act as if S were parallel to. the 
first principal axis. 

Since Y= Q;—0 is a consistent subsidiary condition, we may obtain 
the equations of motion of the remainder of the variables by setting Q,=@3=0 
in the hamiltonian. The parts, H,-+H,, then reduce to the Pauli hamiltonian 
(see paper (A) eqs. (31) and (35)). The term, H, becomes 


where we have written Q, = ho/2. This part of the hamiltonian will have 
the same effect as adding to the Hamilton-Jacoby equation (paper (A), eq. 
(34-b)) a «rest mass» term such as would be obtained by deducing the Pauli 
equation as the non-relativistic limit of the Dirac equation. To bring out 
this property of H, more directly, we may make the canonical transformation 
to the Pauli spinor, A, given in eq. (22), which yields the «rest mass » term, 


, Nh? e 
(33) He sr | (AA) Ax. = 

We conclude then that if the subsidiary condition, Q, = Q, = 0 is satis- 
fied, our set of equations representing the motion of a field of spinning bodies 


ta 
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will reduce essentially to the Pauli equation, which is a linear equation for 
a spinor, Ÿ. Moreover, even if the spin does not point along the principal 
axis, the components Q, and @; will in general turn with some angular ve- 
locity, ©, so that for processes involving characteristic times, 7 >1/w, the 
Pauli equation will still serve as an adequate approximation, since the effects 
of the components Q, and Q, will average out to zero in such processes. But 
in processes so rapid that t= 1/«, the Pauli equation would in such cases 
cease to be a good approximation, and one would have to go back to our 
basically non-linear set of equations for the various spin and angle variables. 
Of course, in the present theory, there is nothing to determine the value of ©. 
We may, however, consistently suppose that 1/m is of the order of the times 
involved in the high energy processes connected with the «creation », « des- 
truction » and «transformation » of so-called «elementary » particles. In this 
case, the linear Pauli equation would be a good approximation for atomic 
processes, but would break down completely with regard to the high energy 
processes described above. Of course, to treat such processes correctly, we 
should need a relativistic theory, and this will require a causal interpretation 
of the Dirac equation. But we can already see that the generalization of the 
Pauli equation adopted here leads in a natural way to a breakdown in con- 
nection with high energy processes of the hypothesis of linear superposition, 
which is one of the basic postulates underlying the present form of the quantum 
theory. When this postulate fails, then the usual interpretation of the quantum 
theory cannot consistently be applied. This is an example of how a causal 
interpretation of the quantum theory permits the consideration of new kinds 
of theories, not permitted if one restricts oneself to the theories that are con- 
sistent with the usual interpretation. 

Finally, we note that in the present theory, the basic meaning of # is 
that 7/2 is the angular momentum per body in the fluid. Thus, in the theory 
proposed here, the statement that # is a universal constant of nature implies 
that all bodies have the same spin angular momentum. This requirement is, as 
we have seen, a subsidiary condition that is consistent with the equations of 
motion that we have adopted. Of course, we have not yet explained why 
the subsidiary condition should be universally satisfied but this may perhaps 
be done later in connection with more extensive theories, such as, for example, 
those concerned with a causal interpretation of the Dirac equation and second 


quantization. 


4. — An Illustration in Terms of Orbits in a Hydrogen Atom. 


We shall now illustrate the causal interpretation of the Pauli equation 
in terms of the example of an electron in a hydrogen atom. This illustration 
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will also permit us to draw some interesting conclusions concerning the meaning 
of the quantization conditions in the hydrodynamic model with spin. 
We begin by considering a stationary state, in which the spinor wave 


funetion «, is given by 


(34) o = R(x) exp [— iHt/h] de si 


Be (x) 
cos 0/2 exp [ip/2] | 


= R(x) exp [iw(x, #)/2] 
à sin 0/2 exp [— ip/2], 


2 js the normalized spinor of eq. (7) and 


where R(x) — V jo |? JE | 


(35) ye, )=—7 pa), with ya) =[yl@, 0]. 
(Note that here y(x, t) refers to the Euler angles of the body which is a function 
of position, while « is the wave-function spinor). 

We see then that in a stationary state the orientation angles 6 and are 
in general functions of position, but are constants in time. However, the 
angle y of rotation about the first principal axis increases linearly with the time, 
when evaluated at a given position, x. 

If, however, we look at a particular moving body, then it will follow some 
orbit in space, in which the various Euler angles and spin directions may 
change with time. But the motion must be periodic, in the sense that after 
the body returns to its original position in space, the Euler angles 0 and 
must return to their original values (plus perhaps some multiple of 2x for 
the Euler angle g) while y must have changed by — 2HAt/h (plus a suitable 
multiple of 27). If these conditions were not satisfied, then we could not have 
a stationary state in which the only change with time of the state of the fluid 
at a given point is at most a rotation of each body around its axis of symmetry. 

To show how this wave function leads to quantization of angular momentum, 


we shall now evaluate the action integral, 7 = Pp ‘dx, where p is the mo- 


mentum of the motion of the center of mass of the body in its periodic orbit. 
Expressing p with the aid of paper (A) eq. (23), we obtain 


h 
pidre 5 (Vp + cos 0 Vo): dx 


Now dx represents the change of x as we follow the moving body in its orbit. 
Let us then introduce the differential, 6x, which represents the change of x 
at a fixed instant of time. Now, from eq. (35), we have Vy = Vy (x). Since 
neither 0, g nor y, are functions of t, we can then replace dx in the action 
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integral by dx, obtaining 


SR h 7 
(36-a) p'dx = = (Vy + cos?Vp): da = d (dy) + cos 0 dp) 
and 
; : xl 
(36-b) D p'Ax— È Cp (dy + così 0@) , 


where on the right hand side of (36-5), the integration is carried out around 
a circuit in space taken at a particular time, but the same circuit as is made 
by the actual orbit of the body. 

We must now find out how y, and q change as we go a around the circuit. 
Obviously we must have Ay) = 21,2 and Ag = 2n.% where n, and n, are 
integers, in general different from each other. But now we shall show that 
n, and n, must satisfy the further conditions that n,+n, and n,— M: shall 
be even numbers. 

To prove this, we must first briefly review some of the topological properties 
of the rotation group (°°). Consider, for example, a series of infinitesimal 
space displacements in our field of body orientations, which add up to give 
some finite displacement. Such a series of space displacements will lead to 
a corresponding series of infinitesimal rotations, which add up to give the 
total rotation needed to bring the body from the orientation that it had at 
the beginning of the series to the one that it has at the end. Now, if the 
series takes the form of a closed circuit, it is evidently necessary that the body 
finally return to its original orientation, so that, as we have already pointed 
out, an integral number of rotations of 27 must have taken place along the 
circuit. But now, let us suppose that we make the circuit smaller and smaller, 
permitting it to shrink continuously down to a certain point. The total ro- 
tation that takes place as we follow the circuit must then decrease continuously 
to zero. Otherwise, we should have a finite rotation connected with an infi- 
nitesimal displacement in space, and this would not be consistent with the 
requirement that the field of body orientations be continuous. 

In the mathematical expression of the above requirement of continuity 
in the field of body orientations, the Buler angles do not provide a convenient 
parametrization, because when 0 = 0, singularities can appear in the para- 


(5) H. WEYL: Theory of Growps and Quantum Mechanics (New York, 1928). See 


p. 180. 
(6) F. D. MURNAGHAM: The Theory of Group Representations (Baltimore, 1938). 


See p. 318. 
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metrization that do not correspond to any real physical singularity in the 
field of orientations. Indeed, when 0 = 0, the only quantity defined by the 
actual rotation of the body is the sum of the angles, y+y, while the difference, 
p— y is completely arbitrary. Thus, if we deform our circuit past a point 
where 0 = 0, it is possible for y and y separately to undergo large changes, 
even when the actual rotation (and therefore along with it the sum, +) 
changes only by infinitesimal amounts. Hence continuous changes in the 
orientations of the bodies need not always be reflected as continuous changes 
in the Euler angles; and for this reason, there is no guarantee that as we shrink 
a circuit in a field of rotations down to a point, the corresponding circuit in 
the space of the Euler angles will also shrink down to a point. As a result, 
the expression of the requirement that the field of body orientations be every- 
where continuous would be rather complicated if done in terms of Euler angles. 

To avoid the kind of difficulty described above, we may parametrize the 
rotation group, not in terms of Euler angles, but rather, in terms of the Cayley- 
Klein parameters, which are (according to paper (A), eq. (8-b)): 


b,= cos0/2 cos(w+ ¢)/2, b, = cos 0/2 sin (y + ¢)/2, 
bs = — sin 6/2 sin (w— @)/2, b, = sin 0/2 cos (y — p)/2. 


For the Cayley-Klein parameters have the advantage that they provide 
everywhere a locally unique and continuous representation of the rotation 
group, in the sense that an infinitesimal rotation always implies a corresponding 
infinitesimal change in these parameters. To demonstrate this property of 
the representation, we note that the unit spinor (which is in a one-to-one 
correspondence with the Cayley-Klein parameters) transforms under the infi- 
nitesimal rotation, % dt, as 05 = i(dt/2)(W-0)8. Thus, an infinitesimal ro- 
tation must always produce a corresponding infinitesimal change in f (and 
therefore in the C.K. parameters). 

It must be remembered, however, that the unique character of the C.K. 
parameter representation applies only to small rotations. Indeed, with regard 
to the properties of the rotation group «in the large », this representation is 
two-valued. For corresponding to any given rotation, there always exist two 
sets of the b,, one of which is the negative of the other. But even the two-valued 
character of the C.K. parameters in the large is significant; for it provides 
a mathematical description of the fact that the rotation-group is doubly- 
connected, a fact that is, as we shall see presently, very important in the de- 
termination of the changes of angle that are permissible in going around a 
circuit. 

To see the meaning of the two-valued character of the representation of 
the rotation group in terms of the C.K. parameters (and therefore in terms 
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of the spinors), let us consider a circuit in space that leads to a continuous 
series of infinitesimal rotations of the body orientations. There is of course 
an ambiguity in the point in the C.K. space at which we must start, because, 
corresponding to the initial orientation of the body, there are two possible 
sets of values of the C.K. parameters. But after we have chosen one of these 
sets, then the path in C.K. space that corresponds to our series of rotations 
is determined uniquely because of the locally continuous character of the 
representation of the retation group in terms of the C.K. parameters. Let us 
now consider a circuit leading to a rotation of 27 about some axis. It is clear 
from the definition of the b, that such a rotation carries us from a given set, 
(b,, bs, b3, ba) of C.K. parameters to the corresponding set (— di, DEA bs, 
— bi). Thus, even though the body has come back to its original orientation, 
the C.K. parameters have not come back to their original values. Indeed, 
geometrically speaking, this rotation carries us only half way around the hyper- 
spherical surface on which we must remain in the four dimensional space of 
the b;. As in the case of a great circle on a sphere in three dimensional space, 
there is no way to continuously deform such a curve to a point while keeping 
the endpoints fixed (7). Asa result, when we make a closed circuit in the field 
of body orientations, if the body rotates through 27 along this circuit, then 
as the circuit is shrunk down to a point, the total rotation that takes place 
along this curve will have to remain equal to 27, no matter how small the 
circuit becomes, so that the field of body orientations cannot be continuous. 
On the other hand, if the circuit carries us through two rotations of 27, we 
are brought back to the original values of the C.K. parameters; and in the 
space of these parameters, the corresponding path goes all the way around 
the hypersphere. Then, as in the case of a complete great circle on a sphere 
in three dimensional space, this curve can be deformed continuously to a 
point, while keeping the endpoints fixed. From this, we conclude that a 
circuit involving two rotations of 27 can be shrunk continuously to one having 
no rotation at all, while a circuit involving a single rotation of 2x cannot. 
More generally, this will happen wherever g-+v and gp y change by even 
multiples of 27 on going around à circuit; or whenever "+ and n, — Nez 
are even integers. Thus, we justify the conditions which we gave earlier, re- 
garding the permissible changes of value that the wave function may suffer 


on going around a circuit. 
We may illustrate the conclusions of the preceeding paragraph in terms 


of an example in which a given circuit carries us through two rotations of 27 


(7) The fixing of the end points corresponds to the fact that we are choosing a 
circuit that passes through a certain point in space. Since the orientation of the body 
is fixed at this point, the C.K. parameters cannot change at the endpoints of the 
circuit, no matter how this circuit is shrunk or otherwise changed. 
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in any specified direction. Now, as the circuit is shrunk, it would be possible 
for the rotation to become equivalent to two separate rotations of 27, each 
carried out in a different direction; for in this case, the initial and final orien- 
tation of the body would still remain unchanged when the circuit was altered. 
Thus, a continuous deformation of the path in the C.K. space would take 
place, with the endpoints of the path held fixed; for the direction of one of the 
rotations of 27 could change continuously relatively to that of the other. 
But when we arrived at a path in which each rotation had an opposite sense, 
to that of the other then the two rotations would cancel each other; and we 
would have no rotation at all. Thus, we have given an example of how a 
circuit involving two rotations of 27 can be deformed continuously into one 
involving no rotation at all. As we have shown, however, this cannot be done 
with a circuit involving a single rotation of 27. 

We see then how the two-valued character of the C.K. parameters (and 
therefore of the spinors) reflects the topological connectivity properties of the 
rotation group. Of course, if we had to deal only with a single body, these 
topological properties would not be relevant, because a single rotation of 27 
would, after all, bring the body to an equivalent orientation in space, so that 
at a given time, no physical difference could exist between a body that had 
suffered only a single rotation of 27, and one that had suffered two such ro- 
tations. On the other hand, because we are dealing with a field of body 
orientations, the relative amounts of rotation suffered by bodies in different 
parts of the field become relevant; in the sense that if the orientation changes 
by an odd number of rotations of 27 along a circuit, then this circuit cannot 
continuously be shrunk to one with no rotation at all, while if it changes by 
an even number of rotations of 27, it can so be shrunk. Thus, the two-valued 
character of the spinors (and of the C.K. parameters) describe physically signi- 
ficant properties of a field of body orientations. In this way, we obtain a 
reason why the basic quantum-mechanical theory of the electron which con- 
cerns itself with such a field of body orientations, should be expressed in terms 
of spinors. 

We now apply the above conditions on the continuity of the angles of 
rotation to a stationary state of a hydrogen atom, which in the usual inter- 
pretation is described as having an orbital angular momentum, J, and a Z 
component of the total angular momentum (spin plus orbital) of X — (1 + 4). 
For this state, the wave function is (8) 


. VIE PE (0 xp [i(k — 3)g' 
3 = Re) | di di (0') exp [i(k —4)q'] | 
V2 +1 VI 5 1 PH (0)  expli& +9] 


(5) See D. Boum: Quantum Theory (New York, 1951); Chap. 17, eg. (ta): 
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where 6’ and g' are the polar angles of the position vector, x, of the center of 
que (k—}) =. . > > 
the body, while P* >(0’) is a suitably normalized associated Legendre poly- 
nomial. 
Comparison with eq. (34) and (35) then indicates that 


(38) eee ig! = ne RUN? 
BAER) 

Since }—1-+ 4, and since g' changes by 27 as we go around a circuit, 
we see that Ay, = (21 + 1)27 and dgp=— 2x. Thus, the topological rela- 
tionship that m+ y and p— y change by even multiples of 27 on going around 
a circuit is satisfied. | 

We are now ready to compute the action integral (75) for this case. To 
do this, we note that the velocity, v—(h 2m)(Vy + cosO Vo), has only a com- 
ponent in the direction of g'. Thus, the body moves in a circle, the plane of 
which is normal to the Z axis. From (77), we see that 0 is constant along 
this circie. Thus we obtain 


208 1 
(39) fea) (i eee " aap È p_(1--cos 9); 


€ 
2 2 


\ 


where n is an integer lying between 1 and — 1. 
The Z component of the orbital angular momentum of the body is 


T 108 0 
(40) Pe = ae =F (h SE ) 


The combined Z component of spin and orbital angular momentum is 


h cos 0 ; 
(41) DE ms hk . 


Thus, we find that the total angular momentum is just equal to that pre- 
dicted in the usual theory. This angular momentum divides itself, however, 
between the spin and orbital angular momentum in accordance with the co- 
latitude angle 0’ of the radius vector of the body relative to the center of the 
atom (which determines 0 through eq. (38)). 

The result (39) for the action variable is analogous to that obtained from 
the Bohr-Sommerfeld quantum condition (of the «classical » quantum theory ) 
in which even the feature of fractional quantum numbers is present, because 
these are often found to give better agreement with experiment than can be 
obtained with integral quantum numbers. There are, however, three important 
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differences from the Bohr-Sommerfeld theory. First, eq. (39) is an exact 
“equation, and not an approximate one. Secondly, in the evaluation of J, the 
momentum, p, must be calculated taking into account the quantum-potential 
and the spin energy, and not just the classical potential as was done in the 
Bohr-Sommerfeld theory. Finally, the action variable, J, for the particle 
motion alone depends in the location of the orbit. Only the sum of this action 
variable and the one corresponding to the spin angular momentum is quantized 
in a way that is independent of the location of the orbit. 

We are now in a position to see how our model accounts for the quantization 
of the angular momentum. To do this, we first note that the action variable 


is proportional to p v-dx which is an integral representing the total vorticity 


of the fluid inside the circuit in question. Thus, we have obtained a simple 
physical interpretation of the action variable. Then, because of the relationship 
v = (h/2m)(Vyp + cos0Vy), the connection (36) is established between this 
vorticity, and the changes of body angle, Ay, and Ag, which occur as we go 
around a circuit. Finally, because we have a continuous field of spinning 
bodies, 4y,-+4g and Av— Ag must be even multiples of 27. Thus, the 
quantization of vorticity and therefore of angular momentum reflects basically 
the requirement that we have a continuous and single-valued field of spinning 
bodies, whose orbital motion is coordinated to their spin motion in the way 
implied by the relationship, v = (h/2m)(Vy + cos 0 Vo). 


5. — On the Process of Measurement of the Spin. 


We shall now indicate in general terms how the process of measurement 
of the spin is to be treated in the model given in this paper. To do this, we 
shall first review briefly a similar treatment of the analogous problem of a 
non-spinning particle (described by the Schrüdinger equation and not by the 
Pauli equation) going around an atomic nucleus with orbital angular mo- 
mentum, #. 

Now a common method of measuring the spin of an atom is to place the 
atom in a non-homogeneous magnetic field. This field then separates atoms 
according to their component of the angular momentum along the direction 
of the field. To show how this process of measurement is treated in terms of 
the causal interpretation of the quantum theory, we shall use a method de- - 
veloped in a previous paper (*). Now the general initial wave function of the 
electron with orbital angular momentum of # is 


(42) P = ax) + Pe) + di P(x), 


() D. Boum: Phys. Rev., 85, 180 (1952). 


— 
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where ¥Y_,, %, W, represent the wave functions of electrons with Z com- 
ponents of the angular momentum of respectively, — #, 0 and h, while a_,, 
M, 4, are the corresponding coefficients (in general complex) for the expansion 
of the wave function, (x). The complete wave function for the combined: 
system consisting of electron plus the atomic nucleus (whose coordinates we 
denote by y) is then 


(43) foly)(a_.P_.(x) + ap¥o(x) + ty Fo (x) ’ 


where f,(y) represents a wave packet describing the fact that the atomic 
nucleus is fairly well-localized in space. 

Now when we apply a magnetic field, $, (the direction of which we take 
to be that of the Z axis) then each of the three parts of the wave function 
described above will begin to oscillate at a different frequency. Moreover, 
if the magnetic field is inhomogeneous, the three parts will begin to separate 
in space. Indeed, after some time, the wave function will be transformed into 


(44) P= fly A)a_, exp [ia ]P_1(2) + fo(Y )Ao EXP Lio] Po(x) + 


lo Vins A)a, exp [ix] (a) ; 


where x_;, %, % are changes of the phase angles which have resulted from the 
effects of the magnetic field, and A is the motion of the nucleus under the 
action of this field. Now eventually A becomes much bigger than the width 
of the wave packet f,(y), so that the three parts of the wave function cease 
to interfere with each other, and obtain a classically distinct separation. When 
this happens, then as has been shown (*), the particle-like inhomogeneity 
associated with the electron must have entered one of the separate packets. 
The probability that it is any one of them is given by 


2, Pili 


(45) SERE doa: 4 Py =| Ge 


Thereafter, the other packets play no role, so that they can be neglected. 

We see then that what really happens in a measurement is that the meas- 
uring apparatus provides a general environment (in this case, the inhomo- 
geneous magnetic field) in which the wave function is transformed from what- 
ever it may have been initially into an eigen-function of the « observable » 
that is being measured (in this case, the Z component of the angular mo- 
mentum). 

Thus, in the case under discussion, the actual orbital angular momentum 
of the electron can vary continuously but only certain stable quantized values 
are possible, which can exist indefinitely without change, in the environment 
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supplied by the inhomogeneous magnetic field. Which of the three possible 
stable values will actually be obtained is not determined in an individual case, 
if we merely specify the initial wave function, (x). To determine, the actual 
result, we should have to specify also the initial location, §(1,), of the particle, 
and the general irregular motions in the Ÿ field which lead to a statistical 
behavior that is described by the probability distribution, P —|W|*? for the 
particles. But as we have seen, the probability of obtaining any particular 
result in a statistical aggregate of cases is determined in terms of the coeffi- 
cients of the initial wave function with the aid of eq. (45). 

In the causal interpretation of the quantum theory, a typical measurement 
process of the kind that can be carried out in connection with the atomic 
level is therefore not really a measurement of the detailed properties of the 
underlying system (which is assumed to consist of the Ÿ field plus the particle 
that moves in it). Rather, it represents a kind of statistical response to certain 
over-all properties of the Ÿ field and of the particle. In this sense, the « ob- 
servables » are rather analogous to the pressure and temperature of macro- 
scopic physics or to the macroscopic variables used in hydrodynamics, since 
these likewise do not describe any detailed properties of the underlying mo- 
lecular motions, but rather general statistical properties of the system as a 
whole. 

Let us now consider how this theory of measurements could be extended 
to the Pauli spin theory. To do this, we should need a theory of the many- 
body problem, corresponding to our theory of the many-body-Schrédinger 
equation (°), because as we have seen, we have had to discuss the interaction 
of the electron with the atomic nucleus in order to develop a theory of measu- 
rements. Such a theory is now being developed in terms of an extension of 
our model to include second quantization; and preliminary results suggest 
that such an extension will be possible, although many of the details remain 
to be worked out. For the present, however, we shall merely state that it 
appears that the theory of measurements can be carried out for the spin theory 
in a way that is essentially the same as what has been done for the Schrô- 
dinger equation without spin. For example, in the case of a spinning electron 
in an atom in an «s» state, there will be two possible basic solutions, one 
corresponding to spin up, and the other to spin down, with the spin measured 
in any desired direction. Hence, the initial wave function can be expressed as 


Pi ay x) a Aes 


Then, as in the case that we have already discussed of a particle of angular 
momentum, #, obeying Schrédinger’s equation, the effect of the inhomogeneous 
magnetic field will be to transform this wave function either into ‘7. or into 
‘P_, with respective probabilities, P.,=|a+|?, and P_=|a_|*, of obtaining 


se 
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either result. Hence, as in the Schrédinger theory, the appearance of quantized 
possible values for the spin in a measurement will be the result of a reaction 
to the magnetic field of the system as a whole, (consisting of fluid with spinning 
bodies, and the particle-like inhomogeneity, described in Sec. 1, which is the 
counter part of the particle appearing in Schrédinger’s equation). Thus, what 
is obtained in what is now called a measurement of the spin will have no 
direct and simple relationship to the spins of the bodies constituting the fluid. 
Indeed, as we have seen, the components of the latter can vary continuously, 
but nevertheless, in a magnetic field, the overall motion of the whole system 
is such that the total angular momentum eventually settles down either to 
h/2 or to —h/2. (We may make here an analogy to certain kinds of classical 
non-linear oscillators, which after being disturbed eventually settle down to 
one of a number of possible stable modes of oscillation). It is clear then that 
the spin as it is now measured should be considered as a higher-level property, 
having a relationship to the assumed spinning bodies that is somewhat ana- 
logous to the relationship between macroscopic variables, such as pressure 
and temperature, and the underlying atomic variables. 


6. — Summary and Conclusions. 


In this paper and in the previous paper (A) we have developed a model 
for the Pauli equation in terms of a fluid composed of spinning bodies, which 
contribute an «intrinsic angular momentum » to the total angular momentum 
of the system. This model has the property that if the bodies are at any time 
all spinning with their angular momenta parallel to their principal axis of 
symmetry, then they will continue to satisfy this condition for all time. On 
the other hand, it is possible for the angular momentum S to have a general 
orientation; and in this case, the component of S normal to the principal axis 
will turn with an angular velocity, ©, that depends on how fast the body 
happens to be spinning and in the torques acting on the body. For processes 
with characteristic times, t > 1/0, the component of S normal to the prin- 
cipal axis will average out to zero, and the Pauli theory will provide a good 
approximation. But for processes in which 7 is of the order of 1/w or less, 
the Pauli equation will no longer apply, and the full general set of non-linear 
equations will be needed. This means that our model already implies the 
possibility of a break down in connection with sufficiently high frequencies, 
and therefore with sufficiently high energies, of the whole general scheme 
connected with the usual interpretation of the quantum theory, which is based 
in an essential way on the assumption that the fundamental equations of the 


theory will always be linear. 
In connection with our discussion of the theory of measurements in Sec. 5, 
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it was seen that the spinning bodies of which our fluid is assumed to be 
constituted are not identical with the spin « observables » of the usual quantum 
theory, but that rather, they constitute a lower level, in terms of which the 
spin «observables » are determined as overall and in general statistical pro- 
perties of the fluid. For example, the characteristic quantized way in which 
the spin angular momentum manifests itself at the atomic level was seen in 
Sec. 4 to follow from conditions of single-valuedness applying to the motion 
of the fluid as a whole, which are such that even though the spin motions of 
the bodies are continuous, the overall motion has certain discrete possible 
stationary values for the angular momentum. 

In sum, then, it may be said that we have, for the case of the Pauli equation, 
explained the quantum theory in terms of the motions of new entities existing 
at a sub quantum-mechanical level. We call the new level «sub quantum- 
mechanical » because the laws of quantum-mechanics do not apply there. 
Rather, the laws of quantum mechanics emerge as overall and statistical 
relationships arising on the basis of the lower level laws, as for example, the 
laws of ordinary hydrodynamics arise on the basis of lower level laws governing 
the atomic motions. 

Naturally, to make an explanation of the quantum theory possible, we 
have had to postulate something, viz, the fluid composed of spinning bodies 
(since without assuming something we can never explain anything). It may 
then be asked what we have gained by making such a postulate. First of all, 
we have gained the possibility of seeing in a rational way how all of the phe- 
nomena of atomic physics could be connected by means of a set of general 
causally determined motions, so that we do not have to regard atomic pheno- 
mena as mysterious processes which take place in a way that could never 
even be conceived of. (In this connection, let us recall that there has existed 
a widespread general impression that to obtain such a rationally understand- 
able explanation of quantum phenomena in general and of the electron spin 
in particular would be impossible.) Secondly, whenever one obtains a rational 
explanation of a wide range of phenomena past experience in science has shown 
that this explanation generally suggests fruitful new avenues of approach to 
problems, which would not even have been suspected if the phenomena had 
not been thus explained, but rather had simply been accepted as things that 
« just happen » for no particular reason whatever. For example, the atomic 
theory, suggested originally by the effort to explain the laws of chemical com- 
bination and the gas laws in terms of the properties of atoms, was eventually 
able to explain many new kinds of phenomena (e.g., Brownian motion, visco- 
sity, gaseous discharges, etc.) and suggested important new directions of re- 
Search (e.g., Rutherford scattering, electron theory of metals, ete.). As for 
the question of wether the assumption of a sub quantum-mechanical level 
will eventually prove to be fruitful in a similar way, this can of course be 
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answered definitively only in the future. Nevertheless, one can already see 
good reasons why this approach may be on the right track, even if, perhaps, 
not correct in all of its details. Thus, the characteristic new phenomena of 
modern high energy physies is the appearance of a whole host of «element- 
ary » particles, which can be «created », «destroyed » and transformed into 
each other. The very fact that these processes of creation, destruction and 
transformation are possible suggests strongly that the so-called elementary 
particles are not really elementary, but rather, that they arise on the basis 
of motions of new kinds of entities that are still more fundamental. Thus, 
what is suggested is a new level, below that of the « elementary » particles, out 
of which these particles arise as some kind of moving structures. 

Now, we have already seen that to explain the quantum theory causally, 
we have already had to postulate a sub-quantum mechanical level, out of 
the motions of which the usual quantum mechanical properties of things arose 
as overall characteristics (e.g., quantization). Now, as long as the basic 
equations governing the system are linear, nothing new can arise in these 
overall characteristics, not already treated in the well-known solutions of the 
Pauli equation. But as we have seen, it is just in connection with sufficiently 
high energy processes that the equations of our model can become non-linear. 
Now, it is well known that non-linear equations have, in general, many modes 
of stable motion. Each of these modes would manifest itself at the atomic 
level in connection with new rules for quantization and for the determination 
of other overall properties of the system, which we would interpret in terms 
of the appearance of a new kind of « particle ». Thus, the way is opened up 
for a treatment of the processes of creation, destruction, and transformation 
of elementary «particles », as well as for a calculation of which kinds of 
« particles » can exist, and of what some of their properties are, since the new 
« particles » could correspond to new modes of overall motion of the underly- 
ing fluid. 

Of course, we do not believe that a model based on an explanation of the 
Pauli equation will really be adequate for the purposes described in the previous 
paragraph because it is not relativistic. A model based on an explanation of 
the Dirac equation (and better still with second quantization) should however 
give a much more accurate treatment than would be possible with the model 
given in this paper. Present work indicates.that models can already be found 
which reproduce most of the features of the Dirac equation and many of those 
of second quantization. The completion of this work would then lay the 
foundation for an attack on the properties of the new level, including those 
connected with the creation, destruction and transformation of elementary 
particles. In any case, it is clear that new directions of investigation could 
thus be opened up, going outside the framework of theories that fit into the 
current general scheme of the quantum theory. 
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